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I.  Summary 

Simulation  methodologies,  algorithms,  and  computer  codes  allowing  molecular  dynamics 
simulations  of  formation,  propagation,  and  breakup  processes  of  nanojets,  generated 
either  through  the  application  of  pressure  or  through  the  action  of  an  electric  field  to 
liquids  containing  solvated  charges,  were  developed.  Particular  emphasize  was  placed  on 
simulations  of  electrically  driven  nanojets  of  liquid  formamide  solutions  with  dissolved 
Nal,  and  on  investigations  of  the  behavior  and  response  of  nanoscale  drops  of  such 
solutions  as  well  as  of  pure  formamide,  when  placed  in  strong  electric  fields.  The 
simulations,  and  coordinated  electrospray  experiments  at  the  the  AFRL,  Hanscom  AFB. 
demonstrated  a  mixed  charge  emission  regime  exhibiting  field-induced  cluster  ion 
evaporation  and  ejection  of  charged  droplets.  The  simulation  average  mass-to-charge 
ratios  and  maximal  surface  fields  of  about  1  V/nm  agree  with  the  experimental  results  and 
with  electrohydrodynamic  theory  of  cone-jets.  The  measured  solvated  ion  distributions 
are  also  correctly  reproduced  by  the  simulations.  Emission  of  charged  particles  is  found 
in  the  simulations  to  occur  predominantly  from  the  tip  of  the  nanojet,  rather  then  from  the 
neck  between  the  Taylor  cone  and  the  nanojet.  The  MD  simulations  revealed  novel  field- 
induced  structural  ordering  processes  and  electro-crystallization  of  pure  formamide 
drops.  The  microscopic  results  obtained  through  the  atomic-scale  simulations  were 
analyzed  and  compared  to  the  predictions  of  a  continuum  formulation. 

A  detailed  description  of  the  molecular  dynamics  simulations  of  nanojet  formation,  field 
-  induced  cluster  ion  evaporation  and  charged  droplet  ejection  is  given  in  section  IV, 
titled: 

NANOJETS,  ELECTROSPRAY,  AND  FIELD-ION  EVAPORA  TION:  MOLECULAR 
DYNAMICS  SIMULATIONS  AND  LABORATORY  EXPERIMENTS. 
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This  section  gives  a  summary  of  some  of  the  main  results  obtained  through  the  molecular 
dynamics  simulations  performed  in  the  Center  for  Computational  Materials  Science  at 
Georgia  Tech,  together  with  a  description  of  coordinated  laboratory  experiments 
performed  at  the  AFRL,  Hanscom  AFB,  under  the  direction  of  Dr.  R.  A.  Dressier. 
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VI.  Nanojets,  Electrospray,  and  Field-Ion  Evaporation: 
Molecular  Dynamics  Simulations  and  Laboratory  Experiments 

This  section  contains  along  with  the  results  of  the  molecular  dynamics  simulations 
performed  by  W.D.  Luedtke  and  U.  Landman  at  the  Center  for  Computational  Materials 
Science  at  Georgia  Tech,  experimental  results  obtained  by  Y.-H.  Chiu,  D.  J.  Levandier, 
and  R.  A.  Dressier,  (Air  Force  Research  Laboratory,  Space  Vehicles  Directorate, 
Hanscom  AFB,  MA  01731-3010),  and  electronic  structure  calculations  performed  by  S. 
Sok  and  M,  S.  Gordon  (  Department  of  Chemistry,  1605  Gilman  Hall  Iowa  State 
Universal,  Ames,  JA  500/ 1-31 J 1). 
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1.  Introduction 


Understanding  the  fundamental  properties  of  conducting  liquids  under  the  influence  of 
applied  strong  electric  fields  is  a  topic  of  continuing  and  growing  interest.  The  associated 
phenomenology  governs  a  number  of  venerable  technological  areas,  most  notably 
electrospray  ionization,1  electrospray  propulsion,2  electrospray  painting,  and 
electrospinning  of  materials.  In  all  of  these  applications,  the  electric  field  causes 
deformation  and  charging  of  the  liquid  surface,  with  consequent  phenomena  including 
generation  of  jets,  charged-droplet  formation.  Coulomb  fission,  and  ion-  field  evaporation 
or  desorption. 

Coulomb  fission  of  a  drops  occurs  when  the  forces  associated  with  Coulombic 
repulsion  of  surface  charges  exceed  the  forces  associated  with  the  liquid  surface  tension. 
The  condition  where  the  two  forces  are  equal  is  referred  to  as  the  Rayleigh  instability 
limit.4  Lord  Rayleigh  showed  that  the  charge  of  a  droplet  at  this  limit  is  given  by 

qR  =(8ar2£-0y£/3)i,  (I) 


where  £q  is  the  permittivity  of  vacuum,  y  is  the  surface  tension  of  the  liquid,  and  d  is  the 
droplet  diameter.  The  electric  field  strength  at  the  droplet  surface  in  the  absence  of  an 
external  field  is  given  by 

<1 


E=  E  ■ 
nead 

and  the  cricitical  field,  Er,  is  obtained  from 
9* 


(2) 


4  = 


ne^d 


Jr. 

e0d 


(3) 


As  discussed  by  Dole  et  al.,5  and  subsequently  by  Iribame  and  Thomson,6  large  charged 
droplets  of  volatile  liquids  will  evaporatively  diminish  in  size  (that  is,  reduce  their 
diameter)  while  maintaining  the  charge,  ultimately  reaching  the  Rayleigh  limit  and 
undergoing  Coulomb  fission.  In  the  model  of  Dole  et  al.,5  this  process  repeats  itself  until 
a  charged  solute  residue  is  formed.  Iribarne  and  Thomson  explored  the  conditions  at 
which  the  electric  field  strength  at  the  surface  of  the  droplet  exceeds  a  threshold  value,  £) 
,  at  which  ion  field  evaporation  occurs.  This  condition  can  be  met  provided  E /  <  £)?(t/) 
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(otherwise  Coulomb  fission  occurs),  and  the  solid  residue  diameter  is  less  than  the  critical 
diameter  at  which  E /  =  E^d),  i.e.,  in  the  case  of  salt  solutions  the  concentration  is  low.  Ei 
derives  from  the  energy  to  overcome  a  barrier  associated  with  the  combination  of  the 
solvation  energy  of  the  evaporated  ion,  and  the  electrostatic  interaction  between  the  ion 
and  the  droplet.  The  estimates  of  Iribame  and  Thomson6  put  this  barrier  at  about  4  A 
from  the  surface. 

Based  on  a  kinetic  model  and  charge  mobility  experiments,  Thomson  and 
Iribame7  concluded  that  the  critical  field,  £/,  was  approximately  I  V/nm.  and  that  field 
evaporation  would  be  expected  for  droplets  with  diameters  <10  nm.  Loscertales  and 
Fernandez  de  la  Mora*  and  Katta  and  coworkers9  used  charge  mobility  experiments  and 
refined  models  to  determine  electric  fields  between  0.7  and  1 .9  V/nm  for  the  evaporation 
of  different  singly-charged  ions.  Fenn  et  al.10  also  proposed  similar  fields  by  applying  a 
simple  electrostatic  model  to  the  interpretation  of  dilute  polyethylene  glycol  solutions 
assuming  that  the  charges  of  small  droplets  are  equally  spread  on  the  droplet  surface,  and 
that  the  charge  distances  on  the  droplet  cannot  be  less  than  that  found  on  the  multiply 
charged  large  field-evaporated  ion.  They  determined  fields  ranging  from  1.46  to  2.66 
V/nm  with  increasing  field  strength  required  for  smaller,  less  charged  oligomer  ions. 

Since  Ei  does  not  directly  depend  on  d,  ion  field  evaporation  should,  therefore, 
also  be  a  phenomenon  observed  from  electrified  jets.  The  detection  and  analysis  of 
nascent  ion  currents  from  jets  requires  that  the  spray  is  in  a  vacuum,  and  consequently, 
such  work  is  limited  to  liquids  with  low  vapor  pressures.  Liquid  metals,  having  very  high 
surface  tensions  and,  therefore  high  resistance  to  Coulomb  fission,  are  known  to  produce 
intense  metal  ion  beams  in  an  electrospray  configuration.11  Jet  formation  does  not  occur 
under  typical  liquid  metal  electrospray  conditions. 

Organic  solutions,  the  subject  of  this  work,  have  significantly  more  flexible 
electrospray  properties.  For  example,  with  such  liquids  one  may  observe  a  distribution  of 
emitted  ions,  even  under  conditions  where  jet  and  charged  droplet  formation  are 
favourable.  Most  of  the  early  work  on  electrospraying  organic  solutions  in  a  vacuum  was 
conducted  with  glycerol.  Huberman  conducted  an  energy  analysis  of  positive  ions 
sprayed  from  a  Nal-glycerol  solution,  and  detected  a  feature  that  was  attributed  to 
metastable  solvated  ions  losing  a  solvent  molecule.  Stimpson  et  al.1'  saw  similar 
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behavior  involving  metastables  in  a  mass  spectrometric  analysis  of  both  positive  and 
negative  ions  emitted  from  glycerol  solutions  of  a  number  of  salts.  A  large  number  of 
solvated  ions  as  well  as  protonated  glycerol  ions  were  detected. 

In  the  experiments  by  Huberman  and  Stimpson  et  al.,  the  applied  electrospray 
extraction  voltage  was  very  high,  exceeding  8  kV.  The  high  resulting  fields  are  known  to 
generate  a  stressed  jet  regime  leading  to  multiple  jets  on  the  rim  of  the  capillary  orifice. N 
In  recent  years,  significant  advances  have  been  made  in  understanding  electrosprays 
operating  in  a  substantially  more  characterized  single  cone-jet  mode,  schematically 
depicted  in  Fig.  I.15'24  Taylor25  showed  that  the  surface  tension  and  electrostatic  forces 
for  a  non-conducting  liquid  can  balance  in  a  way  that  leads  to  formation  of  a  liquid  cone 
with  a  well-defined  half-angle  of  49.3°  at  the  cone  apex.  Such  single  cone  configurations 
are  now  referred  to  as  Taylor  cones.  Above  a  critical  field,  the  cone  becomes  unstable 
and  a  jet  is  formed.  Fernandez  de  la  Mora  has  derived  a  more  general  theory  for  single 
cone-jet  sprays  of  conducting  liquids  predicting  a  broader  range  of  cone  apex  angles 
which  was  also  confirmed  experimentally.’'  The  extraction  voltages  to  operate  in  such  a 
Taylor  cone-jet  mode  are  much  lower,  on  the  order  of  1  to  2  kV,  leading  to  a  significantly 
gentler  break-up  of  the  jet.  In  ground-breaking  work,  Fernandez  de  la  Mora  and 
Loscertales"’  derived  both  theoretically  and  experimentally  the  proportionality  between 
the  parameters  of  a  conducting  liquid  and  the  current  emitted  by  the  Taylor  cone 


l~{yKQ!e)'  , 


(4) 


where  K  and  £  are  the  liquid  conductivity  and  dielectric  constant,  respectively,  and  Q  is 
the  volume  flow  rate.  It  is  noteworthy  that  no  electrostatic  parameters  appear  in  the 
expression.  The  role  of  the  external  field  is  the  generation  of  the  Taylor  cone  and  the 
associated  surface-charge  buildup.  Eq.  (4)  can  be  rearranged  in  terms  of  specific  charge 
of  the  emitted  droplets,  q/m. 


(5) 


where  p  is  the  liquid  density.  Since  droplets  formed  in  jet  fission  are  expected  to  be 


charged  near  the  Rayleigh  limit,  the  specific  charge  is  also  a  measure  of  the  diameter. 


7 


Thus,  Eq.  (5)  states  that  the  droplet  size  emitted  by  a  Taylor  cone  can  be  controlled  by 
parameters  such  as  the  liquid  conductivity  and  the  volume  flow  rate. 

Inspection  of  the  cone-jet  structure  in  Fig.  1  allows  the  identification  of  3  regions, 
the  cone,  the  jet  and  the  intermediate  neck  region.  The  latter  represents  the  transition 
from  a  flow  region  in  the  Taylor  cone,  in  which  the  electrostatic  nonequilibrium 
associated  with  the  narrowing  cone  diameter  and  increasing  flow  speed  leads  to  ohmic 
currents  within  the  liquid,  to  a  convective  flow  region  in  the  jet  that  is  in  electrostatic 
equilibrium,  and  where  the  charge  density  is  lowered  as  the  jet  speed  increases.  Thus,  the 
region  of  highest  surface  charge  and  normal  electric  field  is  in  the  intermediate  neck 
region.  Gamero-Castano  and  Fernandez  de  la  Mora22  derived  expressions  for  the 
maximum  normal  electric  field  and  radius  of  curvature  or  jet  radius  in  this  region: 

1  (  K'V 


E  -  — 

Co  '{Q) 


(6) 


4* -KG'*]* 


(7) 


fhe  above  authors  were  also  the  first  to  observe  field  evaporation  from  a  jet  using  Nal 
formamide  solutions.  Volume  flow  rates  of  less  than  100  pl/s  and  conductivities 
exceeding  1  S/m  were  required,  which  lead  to  electric  fields  near  1  V/nm,  fully  consistent 
with  earlier  predictions  of  field  evaporation  onsets.  Conditions  were  achieved  at  which 
the  ion  current  exceeded  the  charged  droplet  current.  At  these  conditions,  the  curvature 
of  the  transition  region  and  the  diameter  of  the  jet  are  predicted  to  be  in  the  vicinity  of  10 
nm.  A  cone-jet  operated  at  these  conditions  can  thus  be  classified  as  a  nanojet. 

Subsequently,  ion  field  evaporation  has  been  commonly  observed  with  a  third 
class  of  low-vapor  pressure  liquids,  namely,  ionic  liquids.26'28  Particularly  low  flow  rates 
can  be  obtained  by  replacing  the  capillary  with  a  needle  etched  to  optimize  wetting 
properties.2^  Pure  ion  emission,  signifying  a  vanishing  jet,  has  now  been  observed  for  the 
ionic  liquids  [EmimJtBFj]28  and  [Emim][Im]  (Emim:  1-3-ethyl-methyl-immidazolium; 
lm:  I-ethyl-3-methylimidazolium  bis{trifluoromethylsulfonIy}imide),30  whereby  the 

former  has  been  successfully  sprayed  using  a  mirco-fabricated  needle  array.31  The 
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present  paper,  with  a  particular  objective  to  learn  about  the  liquid  interface  subjected  to 
high  electric  fields,  however,  will  focus  on  formamide  salt  solutions  because  of  their 
more  interesting  variety  of  field  evaporated  ions,  and  because  they  are  more  tractable 
computationally.  As  demonstrated  by  the  electrospray  ionization  community,  the 
methodology  is  considered  a  so  ft- ionization  approach.  Thus,  ion  field  evaporation  at 
fields  near  the  onset  should  be  a  gentle  process  producing  ions  that  could  reflect  the 
micro  properties  of  the  solution  at  the  interface,  as  was  assumed  by  Fenn  and  coworkers 
in  their  analysis.10 

Recently,  Chiu  et  al.  conducted  mass  spectrometric  analysis  of  positive  ions 
emitted  from  a  concentrated  (28.4%  weight)  Nal-formamide  solution  operated  in  a  Taylor 
cone-jet,  field-ion  evaporation  mode.  The  measurements  identified  primarily  solvated 
ions,  Na+(HCONH2)n,  with  a  preferred  solvation  number  of  n  =  4  which  is  the  largest 
solvated  ion  predicted  to  be  thermally  stable  at  the  room  temperature.  Thus  the 
measurements  appear  to  confirm  the  notion  of  a  thermal  evaporation  process. 
Interestingly,  however,  the  solvated  ion  energy  distributions  were  found  to  be  near  or 
below  the  energies  of  the  charged  droplets,  contrasting  the  electrohydrodynamic 
prediction  that  field  ion  evaporation  would  occur  from  the  neck  region  of  the  jet  which  is 
expected  to  have  a  higher  potential  than  the  tip  of  the  jet  because  of  ohmic  losses  in  the 
cone-jet  transition  region.  Similar  observations  were  made  by  Gamero-Castano,24  for  the 
same  solution.  However,  with  measurements  on  the  ionic  liquid  [Emim][Int]  in  a  mixed 

*  -5 

ion/droplet  regime,  Chiu  et  al.'  ’  could  demonstrate  that  field-cvaporatcd  ions  clearly 
are  emitted  at  the  cone-jet  neck,  as  evidenced  by  their  energy  and  angular  distribution, 
consistent  with  the  theory  advanced  by  Gamero-Castano  and  Fernandez  de  la  Mora.22 

Recent  experimental  and  theoretical  advances'14'26  have  shed  new  light  on  the 
interfacial  properties  of  aqueous  salt  solutions  (under  conditions  free  of  external  electric 
fields).  Counter  to  the  long-held  view  that  atomic  ions  are  repelled  by  the  air/water 
interface  into  the  bulk.  MD  simulations  of  alkali-halide  solutions  have  demonstrated  that 
more  polarizable  anions  reside  preferentially  at  the  liquid-to  vacuum  interface.'4  27-42 
These  results  have  been  confirmed  by  X-ray  photoelectron  spectroscopy  experiments,4' 
as  well  as  sum-frequency29, 40,44' 45  and  second-harmonic  generation  experiments26  4649  of 
aqueous  halide  solutions  at  ambient  pressure.  While  these  studies  examine  the 
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microscopic  behavior  of  static  solution  interfaces,  Vertes  and  co workers50  have 
conducted  MD  simulations  of  an  acidic  positive  charged  water  nanodroplet,  with  the 
amount  of  excess  charge  taken  near  the  Rayleigh  limit.  In  these  simulations  field*ion 
evaporation  of  hydronium  ions  was  observed;  we  note  here  that  in  these  simulations 
negative  ions  were  not  included,  and  this  may  influence  the  properties  of  the  simulated 
system,  particularly  at  reduced  dimensions. 

We  conclude  our  brief  review  of  recent  work  with  the  MD  simulations,  performed 
by  Moseler  and  Landman51,  of  liquid  Hows  at  nanoscale  dimensions,  and  the  formulation 
of  a  stochastic  hydrodynamic  treatment.  In  these  studies  simulations  were  performed  of 
liquid  propane  flowing  under  applied  back-pressure  through  a  cylindrical  tube,  and 
emanating  through  a  nanoscale  nozzle  into  a  vacuum;  recently  Kang  and  Landman  have 
extended  these  simulations  to  injection  of  liquid  jets  into  a  surrounding  fluid  (W.  Kang 
and  U.  Landman,  Phys.  Rev.  Lett.  98,  064504-1  -  64504-4  (2007)).  The  properties  of 
the  resulting  jets  -  in  particular  the  characteristics  of  the  breakup  of  the  liquid  jet  into 
droplets  -  were  found  to  differ  from  continuum  hydrodynamics  predictions.  It  was 
discovered  that  a  newly  derived  hydrodynamic  formulation  that  includes  intrinsic  size- 
dependent  fluctuations,  provides  a  faithfully  accurate  description,  in  agreement  with  the 
atomistic  simulations.  This  work  lays  the  groundwork  for  an  investigation  of  electrospray 
jets,  where  the  cone-jet  transition  region  can  be  regarded  as  a  virtual,  nanoscale  nozzle. 

In  this  paper  we  present  a  systematic  study  of  field  ion  evaporation  from  salt 
solutions  using  a  combination  of  mass  spectrometry  and  MD  simualtions.  To  gain 
insights  into  the  microscopic  processes  of  nanojet  formation  and  charge  emission  from 
dielectric  liquids  under  the  influence  of  strong  electric  fields  we  chose  to  investigate 
theoretically  the  behavior  and  response  of  liquid  drops  (made  of  pure  formamide  or  Nal  s 
formamide  solutions),  since  studies  on  such  system  frees  us  from  complicating  factors, 
including  considerations  of  effects  due  to  flow  rates,  contact  angles,  and  the  geometry  of 
the  capillary  tube  ,or  needle,  used  to  generate  the  slender  liquid  flowing  system.  The  MD 
simulations  were  conducted  on  a  10  am  diameter,  uncharged  (equal  amount  of  positive 
and  negative  charges)  Nal -formamide  droplets  subjected  to  external  electric  fields  with 
strengths  on  the  order  of  1  V/nm.  The  size  of  the  drop  was  selected  such  that  the 
corresponding  curvature  of  the  droplet  (1/R)  would  be  comparable  to  the  estimated 
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minimum  curvature  expected  to  be  encountered  on  a  liquid  jet  (see  Eq.  (7j),  and  thus  the 
physics  of  jet  formation,  ion  field  evaporation,  and  charge  emission  from  such  drops  can 
be  taken  as  representative  of  the  processes  occurring  in  electric  field-driven  cone-jet 
flow  systems.  Indeed,  experimental  results  of  field-induced  dual-jet  formation  on 
droplets,  either  fl owing  through  a  region  of  high  electric  field  or  levitated,  have  been 
reported^'55;  these  studies  relied  on  visual/optical  examination  and  were  conducted  on 
droplets  with  diameters  exceeding  50  pm  exhibiting  formation  of  jets.  Under  these 
conditions  field-ion  evaporation  from  the  main  drop  is  not  expected  to  occur.  We  remark 
here  that  Grimm  and  Beauchamp  conducted  mass  spectrometric  measurements  of  field 
evaporated  ions  from  droplets  of  methanol  solution  that  are  ejected  from  larger  (170  pm 
diameter)  drops.^  However,  the  high  external  fields  that  would  be  required  to  induce 
direct  field  evaporation  from  10  nm  diameter  droplets  (which,  as  noted  above,  is  the  size 
range  of  interest  to  us  here)  were  outside  the  range  of  these  earlier  experiments. 

The  mass  spectrometric  measurements  presented  in  this  paper  were  conducted  on 
ions  that  are  field  evaporated  from  a  cone-jet  of  solutions  (Nal  dissolved  in  formamide) 
with  similar  concentrations  to  those  examined  in  the  simulations.  Quantum  chemical 
calculations  were  carried  out  to  determine  the  energy-optimal  structures  and 
thermochemistry  of  solvated  ion  clusters  with  sizes  equals  to  those  observed 
experimentally.  The  thermochemistry  and  associated  normal  mode  frequencies  are  used 
to  assess  the  thermal  stabilities  and  lifetimes  of  the  ions.  The  solvated  ion  cluster 
lifetimes  are  of  particular  importance  for  our  analysis  because  they  are  used  in 
extrapolating  and  comparing  the  simulated  results  concerning  mass  distributions  of 
emitted  solvated  ions  (that  are  simulated  for  relatively  short  flight  times)  to  the 
experimental,  mass  spectrometric  data,  obtained  for  much  longer  flight  times.  Both 
positive  and  negative  solvated  ion  distributions,  obtained  experimentally  from  respective 
positive  and  negative  extraction  fields,  and  from  calculated  evaporation  statistics  at 
opposite  ends  of  the  simulated  droplet,  are  presented  and  analyzed. 

2.  Experimental 

The  mass  spectrometric  experiment  has  been  described  in  detail  previously,28  and  only 
the  features  relevant  for  this  article  will  be  discussed.  A  schematic  of  the  apparatus  is 
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shown  in  Fig.  2.  The  electrospray  source  consists  of  a  tapered  5  pm  inner  diameter  fused 
silica  capillary  tip  (New  Objective,  Inc.)  that  is  fed  by  a  50  pm  inner  diameter  transport 
capillary  via  a  capillary  union  (Upchurch  Scientific).  The  solution  How  rate  is  controlled 
by  pressurizing  the  liquid  reservoir  supplying  the  transport  capillary  with  dry  nitrogen. 
The  flow  rate  is  determined  from  the  velocity  of  bubbles  introduced  into  the  liquid  near 
the  transport  capillary  entrance.  The  capillary  tip  is  positioned  within  1  mm  of  an 
extractor  electrode  with  an  orifice  of  6  mm  diameter.  The  tip  and  extraction  optics  are  all 
in  a  vacuum  (<10'f>  Torr).  Typical  extractor  voltages  used  for  the  28.4  wt%  Nal- 
formamide  solution  were  2  kV.  The  capillary  bias  voltage  was  applied  to  the  capillary 
through  a  metal  layer  that  coated  both  the  inner  and  outer  surface  of  the  capillary  at  the 
entrance  end.  Typical  bias  potentials  were  450  V,  while  the  extractor  was  biased  at  - 
1,550  V  for  positive  polarity,  and  the  reverse  for  negative  polarity.  These  voltages  were 
chosen  to  minimize  the  ion  energy  in  the  mass  filter,  while  maintaining  reasonable  ion 
transmission  through  the  apparatus. 

The  charged  particle  beam  extracted  from  the  capillary  lip  can  either  be  focused 
on  a  near-field  target  consisting  of  an  interchangeable  Faraday  cup,  or  a  quartz  crystal 
microbalance  (QCM;  XTM/2,  Inficon),  or  towards  the  entrance  of  a  quadrupole  mass 
filter.  The  near- field  target  is  used  for  general  characterization  of  the  charged-particle 
beam  through  total  current  and  mass  flow  measurements.  The  ratio  between  mass  flow 
and  current  provides  the  average  m!q  value,  which,  in  conjunction  with  ion  mass 
distribution  measurements,  provides  a  measure  of  the  fraction  of  charged  droplets  in  the 
beam.  The  near-tield  target  is  translated  away  from  the  beam  axis  for  mass  spectral 
measurements.  Given  the  possibility  of  broad  ion  energy  distributions,  the  quadrupole 
tilter  is  the  mass  spectrometer  of  choice  since  the  quadrupole  mass  setting  is  independent 
of  energy.  However,  the  mass  resolution  declines  with  ion  energy.  Following 
transmission  through  the  quadrupole,  the  ions  pass  through  a  retarding  potential  analyzer 
prior  to  being  detected  using  a  channel  electron  multiplier.  The  retarding  potential 
analyzer  is  used  to  measure  the  energy  distribution  of  mass  selected  ions. 

Mass  calibration,  the  mass  transmission  function  and  detector  response  are 
obtained  from  mass  spectra  of  a  standard  calibrant  compound,  FC-43,  with  a  known 
fragmentation  pattern  at  an  electron  energy  of  70  eV,  The  spectra  were  obtained  after 
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replacing  the  electrospray  source  with  an  electron-impact  ionization  source.  The 
intensities  of  fragmented  ions  of  FC-43  were  measured  for  all  investigated  quadrupole 
transmission  energies,  from  which  the  mass-dependent  discrimination  correction  factors 
were  determined. 

The  concentrated  28.4  wt%  and  1/4  volumetrically  diluted  Nal-formamide 
solutions  used  in  this  study  have  salt  to  formamide  mole  ratios  of  1:8.4,  and  1:36.9, 
respectively.  Similarly,  the  salt/ formamide  mole  ratios  determined  for  KI/Formamide 
solutions  in  this  study  are  1:40.9  and  1:83.0,  respectively.  The  measured  conductivities 
for  the  concentrated  and  dilute  Nal-formamide  solutions  are  1.64  S/m  and  1.26  S/m, 
respectively;  conductivities  of  1 .32  S/m  and  0.77  S/m  are  found  for  the  concentrated  and 
dilute  KI/Formamide  solution. 


3.  Computational  methods 

3.1  Molecular  dynamics  simulations 

To  enable  large-scale  atomistic  simulations  of  the  thermodynamic  and 
electrohydrodynamic  properties  of  complex  fluids  (with  large  partial  and  ionic  charges) 
under  the  influence  of  strong  electric  fields  a  computer  code  with  parallel  architecture 
was  developed.  To  model  these  systems  faithfully,  it  is  important  to  treat  the  long-range 
coulombic  interactions  accurately  and  efficiently.  In  our  initial  studies  of  pure 
formamide  droplets,  long-range  Coulombic  interatomic  potentials  were  treated  with  finite 
interaction  cut-offs  ranging  from  I  to  4  nm.  The  use  of  Unite  range  cutoffs  inhibited  the 
expected  droplet  elongation  as  predicted  by  continuum  theory5*’  of  dielectric  drops. 
Coulombic  interactions  had  to  be  treated  much  more  accurately  (namely  no  cutoff)  to 
bring  the  simulation  results  into  qualitative  agreement  w-ith  existing  continuum  theory. 
Obviously,  the  addition  of  ionic  species  further  increases  the  necessity  of  proper 
implementation  of  long-range  Coulomb  interactions.  Indeed,  the  more  accurate  treatment 
of  dipole-dipole  and  dipole-field  cooperative  interactions  then  allowed  droplets  to 
respond  correctly  to  external  fields  including  field-induced  polarization,  elongation,  and 
instabilities  resulting  in  emission  of  charged  cluster  ions  and  electrified  jets. 
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The  long-range  forces  were  modeled  in  these  simulations  through  the  use  of  the 
fast  multipole  method  (FMM)''’7  specifically  formulated  and  programmed  for  a  parallel 
architecture  and  adapted  to  the  present  problem.  The  FMM  method  enables  the 
modeling  of  large  systems  with  long-range  interactions  by  using  a  “tree-structure” 
subdivision  of  cells  that  allows  a  charge  to  interact  with  a  small  set  of  multipoles 
representing  entire,  increasingly  larger,  volumes  of  remote  space,  rather  than  with  all  of 
the  individual  charges  contained  within  the  cells.  Important  features  of  the  program 
developed  for  this  study  are: 

1)  The  ability  to  replicate  basic  root-cells  and  “stack”  them  along  a  particular  axis  to 
create  a  long  narrow  tree-structure.  This  is  particularly  useful  for  efficient  simulations  of 
systems  such  as  highly  elongated  droplets  and  long  liquid  jets. 

2)  A  large  degree  of  flexibility  in  the  choice  of  boundary  conditions  (bc’s).  The 
computer  program  allows  the  implementation  of  periodic  boundary  conditions  (pbe’s)  in 
one,  two  or  three  spatial  directions,  or  the  execution  of  “cluster  simulations”  where  pbe’s 
are  not  imposed  altogether.  In  simulations  using  pbe’s  in  less  then  3  directions,  the  non¬ 
periodic  directions  can  be  reflecting  (closed  systems)  or  absorbing  (open  system). 
Absorbing  bc’s  (any  atoms  that  cross  the  system  boundary  are  removed  from  the 
simulation)  are  necessary  in  simulations  of  fragmenting  droplets  or  liquid  jets  where  there 
are  large  fluxes  of  atoms  crossing  the  finite  boundaries  (keeping  these  atoms  in  the  root¬ 
cell  would  cause  artificial  effects). 

3)  A  modified  FMM  algorithm  so  that,  when  using  pbe’s,  the  tree-structure  of  the  root¬ 
cell,  along  with  all  its  computed  multipoles,  can  be  replicated  a  specified  number  of  times 
as  neighboring  image  cells.  When  this  is  done,  atoms  in  the  centra!  computational  cell  see 
the  correct  multipole  structure  as  they  interact  with  neighboring  periodic  image  cells  out 
to  any  desired  distance. 

4)  The  ability  to  perform  constant  pressure  simulations.  The  spatial  dimensions  of  the 
root-cell  can  be  allowed  to  vary,  while  using  three-dimensional  pbe’s  and  the  FMM, 
according  to  a  modified  constant-pressure  algorithm  due  to  Berendsenf*  Bulk  samples 
of  both  crystals  and  dielectric  fluids  can,  therefore,  be  prepared  with  solvated  ionic  salts 
at  any  pressure  or  temperature  of  interest. 
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In  the  droplet  simulations  of  this  work,  the  basic  root-cell  of  the  FMM  code  has  8 
tree  levels  with  the  smallest  leaf  cell  chosen  to  be  1 .2  nm  on  each  side.  This  gives  a  root- 
cel!  of  2(8'"  *  1.2  nm  =153.6  nm  dimensions  on  each  side.  For  the  studies  of  droplet 
response  to  a  constant  external  electric  field,  this  root-cell  is  replicated  3  times  along  the 
z-axis  (the  electric  field  direction)  to  give  a  total  tree-structure  of  153.6  nm  x  153.6  nm  x 
460.8  nm  in  the  center  of  which  we  place  the  10  nm  diameter  droplet.  The  very  large 
system  size  allows  enough  room  for  the  parent  droplet  to  move  and  elongate  freely  with 
sufficient  time  for  charged  cluster  emission.  Once  the  simulation  begins,  no  form  of 
temperature  control  or  control  of  the  droplets  center-of-mass  motion  is  performed  to 
minimize  the  risk  of  introducing  artificial  effects. 

Given  the  size  of  the  present  droplets,  it  is  important  to  utilize  an  efficient,  simple 
representation  of  the  inter-atomic  forces.  Formamide  is  treated  as  a  solid  body  with 
quaternion  dynamics,58  implemented  using  a  mid-step  implicit  leap-frog  algorithm'9  that 
has  been  shown  to  be  an  extremely  stable  integration  scheme  with  minimal  energy  drifts 
for  long  simulation  periods.  The  geometry  of  the  formamide  molecule  is  taken  from 
high-resolution  X-ray  studies  of  formamide  crystals,60  The  AMBER  force  field 
parameters61  were  used  to  describe  the  intermolecular  van  der  Waals  interactions  between 
formamide  molecules,  the  interactions  between  the  formamide  molecules  and  the  atomic 
ions,  as  well  as  the  ion-ion  interactions.  CHELP-BOW  partial  atomic  charges62  were 
used  for  the  overall  description  of  the  electrostatic  potential  of  the  formamide  molecules. 
While  long-range  Couloumbic  interactions  between  formamide  molecules  were  not 
truncated,  the  weak  van  der  Waals  interactions  were  truncated  on  a  group  (molecular) 
basis  through  the  use  of  a  smooth  switching  function  described  by  Perera  et  al.63 

3.2.  Quantum  chemical  methods 
3.2.1  Solvated  positive  ions 

We  have  conducted  density  functional  theory  (DFT)  calculations  of  the  observed  positive 
solvated-ion  structures  to  estimate  their  thermal  stability  and  the  dissociation  energies 
associated  with  loss  of  1  formamide  solvent  molecule.  Calculations  were  performed  at 
the  B3LYP/6-31G(d)  level  using  the  Gaussian  98/Gaussian  03  quantum  chemistry 
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package/’4  Geometries  were  optimized,  and  normal-mode  frequencies  for  zero-point 
energy  determinations  were  calculated  for  solvated  ions  of  fNa+(HCONHi)n]  and 
[K+(HCONH2)„],  Where  n=l-6. 

3.2.2  Solvated  negative  ions 

Initial  geometry  optimizations  of  [I_(NH2CHO)  n],  n=l,  5  were  performed  using 
an  augmented  correlation  consistent  double  zeta  basis  set  for  second  row  atom,65  6-3 1G 
basis  set  for  hydrogen,66  Radom  ECP  basis  set  for  iodine,67  and  the  B3LYP68'70 
functional.  These  geometries  were  then  refined  with  Molter-Plessel  second-order 
perturbation  (MP2)  theory71'73  and  the  6-31  lG+(d,p)74  basis  set  using  the  Radom  6- 
31!+G(df)  basis  set  for  iodine.75  Gradient  convergence  tolerance  for  geometry 
optimizations  was  set  to  10'6  Hartree/Bohr.  Stationary  points  were  characterized  by  a 
positive  definite  Hessian  (matrix  of  energy  second  derivatives).  Harmonic  zero-point 
energy  and  temperature  corrections  (to  298. 15K)  were  obtained  using  the  MP2  Hessians. 
All  calculations  were  performed  using  the  June  27,  2005  (R3)  version  of  GAMESS 
(General  Atomic  and  Molecular  Electronic  Structure  System).76  77 

3.3  Unimolecular  ion  dissociation  rate  calculations 

The  distributions  of  the  measured  ion  solvation  numbers  (that  is,  the  distribution  function 
the  number  of  solvating  foramide  molecules  in  a  solvated-ion  clusters)  deduced  from  the 
measurements,  are  compared  with  predictions  based  on  the  assumption  that  the  ions  are 
evaporated  from  a  liquid  at  or  near  room  temperature.  As  the  solvation  number  increases 
the  energy  for  loss  of  a  single  solvent  molecule  reduces,  while  the  internal  energy  of  the 
solvated  ion  increases.  Above  a  certain  solvation  number,  a  fraction  of  the  internal 
energy  exceeds  the  dissociation  energy  of  the  solvated  ion,  making  the  respective  ions 
thermally  unstable.  Using  the  normal  mode  frequencies  from  the  ab-initio  calculations, 
we  determine  the  thermal  internal  energy  distributions  of  the  observed  solvated  ions  from 
the  relation 
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where  Frol(Erol),  and  Fv, *(£«*)  are  the  normalized  rotational  and  vibrational  Maxwell- 
Boltzmann  energy  distribution  functions  for  a  temperature,  7\  and 

Emt  =  Em,  +  Evib.  (9) 

For  those  ions  with  internal  energies  exceeding  the  binding  energy,  we  determine  the 
unimolecular  dissociation  rate  using  a  RRKM  approach,  employing  following  the 
methodology  t  developed  by  Hase  and  coworkers.78  Due  to  the  large  dipole  of 
formamide  (3.7  D)  and  the  resulting  strong  attractive  interaction  between  the  solvated  ion 
and  the  solvating  ligands,  there  is  no  centrifugal  barrier  to  dissociation  and  a  very  loose 
transition  state  is  assumed.  We  assume  conditions  near  the  phase  space  limit,  and  the 
transitional  vibrational  modes  are  treated  as  rotors,  as  prescribed  by  Rodgers  et  al.  All 
systems  considered  include  a  single  solvation  shell  and,  given  the  shallow  global  minima, 
and  near  symmetry  of  the  normal  mode  frequencies,  all  ligands  are  treated  equivalently. 

All  rotors  are  treated  as  active  in  both  the  reactant  and  transition  state,  vibrational  states 
are  treated  as  harmonic  oscillators,  and  the  sum  and  density  of  states  are  determined 
using  the  Beyer-Swinehart  direct  count  algorithm.80  Adiabatic  treatment  of  both  reactant 
and  transition  state  rotors  does  not  markedly  change  the  results. 

The  survival  fraction,  P(nJ),  of  ions  of  a  particular  solvation  number,  n,  at  a 
temperature,  T,  is  then  determined  from  the  expression 

P(n>T)  =  P,(n,T)+  jpjn.T.E^e^’dE^ , 

where  Ps(n,T)  is  the  fraction  of  stable  ions  with  internal  energies  below  the  dissociation 
limit,  Eo\  Pm(n,T,E,„,)dEinl  is  the  fraction  of  ions  with  energies  in  an  interval  Ew,  +  dEm,\ 
k(Emi)  is  the  unimolecular  dissociation  rate  constant,  and  r  is  the  average  Bight  time  of 
the  respective  ions  in  the  mass  spectrometric  apparatus  described  in  Section  2.  The 
experimental  flight  times  are  on  the  order  of  25  x  1 0‘6  s. 


(10) 
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4.  Results 


4.1  Mass  spectrometric  measurements 

Single  cone -jet  sprays  were  generated  for  the  concentrated  and  dilute  Nal  formamide 
solutions.  The  operational  parameters  of  the  experiments  are  summarized  in  Tab,  1.  For 
the  28.4%  wt.  (1:8.4  mole  ratio)  concentrated  solution,  the  extraction  voltage  was  2,000 
V  for  both  polarities,  w'hile  a  slightly  higher  voltage  was  necessary  for  stable  operation  of 
the  dilute  (1:36.9  mole  ratio)  solution.  The  measured  charge  emission  currents 
correspond  to  /(total).  The  ratio  between  the  Faraday  Cup  current,  /(Faraday  Cup),  to  the 
emission  current  provides  a  measure  of  the  beam  divergence.  Mass  llows  registered  by 
the  QCM  detector  are  also  listed.  The  average  charged  particle  mass-to-charge  ratio,  m/q, 
is  then  derived  from  the  ratio  of  mass  flow  to  Faraday  Cup  current. 

Comparison  between  the  total  emission  current  and  that  recorded  on  the  Faraday 
cup  that  collects  a  solid  angle  close  to  that  of  the  QCM  allows  an  assessment  of  the 
electrospray  divergence.  The  measured  divergence  markedly  higher  measured 
divergence  for  the  concentrated  solution  can  be  attributed  to  the  lower  acceleration 
voltage  and  the  higher  current  density,  as  well  as  possible  modification  of  the  emission 
physics  at  the  cone-jet.  This  will  be  further  elaborated  on  in  the  discussion  section. 

The  large  m/q  values  suggest  that  for  both  solution  concentrations  the  cone-jet 
emits  in  a  mixed  droplet/ion  regime,  where  the  main  fraction  of  the  mass  flow  is  carried 
by  the  droplets,  while  a  substantial  fraction  of  the  current  comes  from  ions  generated  via 
the  field  ion  evaporation  mechanism.  This  is  verified  by  a  calculation  of  the  maximum 
electric  field  and  the  corresponding  characteristic  radius  of  curvature  according  to  Eq.  (6) 
and  (7).  The  flow  rates,  Q ,  were  obtained  from  the  mass  flow  assuming  a  density  of  1 . 1 3 
g  cm'1.  The  surface  tension  of  pure  formamide,  0.0575  N  m‘\  was  assumed  for  the 
calculation  of  Emax.  Note  that  Eqs.  (6)  and  (7)  are  based  on  the  total  flow  rate  so  that  the 
derived  values  based  on  the  mass  flow  recorded  by  the  QCM  must  be  regarded  as  upper 
( Emox )  and  lower  {Rmax}  limits.  Depending  on  the  charge  emission  dynamics,  it  is  possible 
that  the  angular  m/q  distributions  exhibit  a  certain  degree  of  non-uniformity.  The 
determined  representative  maximum  electric  field  strength  ranges  from  1.42  and  1.7 
V/nm,  with  higher  fields  corresponding  to  the  concentrated  solution,  while  the  jet 
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dimensions  are  associated  with  radii  in  the  range  between  2.3  and  3.29  nm.  with  the 
larger  radii  determined  for  the  dilute  solution. 

The  negative  total  emission  currents  tend  to  be  almost  twice  as  high  as  those 
observed  for  positive  charge  emission.  Two  significant  differences  are  observed  for  the 
dilute  solution  in  comparison  with  the  concentrated  one:  (i)  for  the  former  the  total 
emission  current  is  lower,  while  the  fractional  current  registered  on  the  Faraday  Cup  is 
higher,  signifying  a  less  divergent  spray,  and  (ii)  the  negative  ion  average  mass-to-charge 
ratio  is  higher  for  the  more  dilute  solution,  while  the  positive  ion  mass-to-charge  ratio 
docs  not  change  significantly  upon  dilution. 

Figure  3  shows  positive  and  negative  ion  mass  spectra  recorded  for  the 
concentrated,  28.4%  wt.  Nal-formamide  solution.  The  spectra  have  been  corrected  with 
respect  to  the  mass  spectral  transmission  function.  The  low  mass  resolution  can  be 
attributed  to  the  high  transmission  energies  of  550  eV.  The  most  abundant  positive  ion 
mass  can  be  primarily  attributed  to  the  solvated  sodium  ions,  Na+(HCONH2)n.  whereby 
the  n  =  4  solvated  ion  is  the  most  abundant  positive  ion  species.  Minor  contributions  to 
the  mass  spectrum,  attributed  to  the  formation  of  Na+(NaI)(HCONH2)m,  are  also 
observed.  The  individual  abundances  are  determined  by  fitting  to  the  observed  peaks 
Gaussian  mass  distributions  given  by  the  experimental  resolution  at  the  expected  peak 
positions.  The  derived  normalized  abundances  are  listed  in  Tab.  2.  Also  shown  is  the 
Na'  solvated  ion  distribution  determined  for  the  1/4  volumetrically  diluted  solution.  The 
abundances  are  essentially  identical  to  those  of  the  concentrated  solution  except  for 
higher  abundances  at  the  highest  solvation  numbers.  Interestingly,  at  the  lower 
concentation  complex  Na+(NaI)(HCONH2)m  ions  could  not  be  detected  with  sufficient 
accuracy  to  allow  a  definite  identification. 

When  applying  an  extraction  field  with  a  negative  polarity,  the  main  observed 
ions  are  solvated  iodine  anions,  T(HCONH2)n>  with  a  distribution  peaking  at  n  =  2. 
Minor  intensities  of  !'(NaI)(HCONH2)m  are  also  observed.  The  abundances  are  listed  in 
Tab.  3.  Contrary  to  the  Na+  positive  ions,  dilution  has  a  marked  effect  on  the  observed 
solvated  V  mass  spectrum,  where  broadening  of  the  distribution  is  observed  with  a  peak 
at  n  =  5.  Similar  to  the  positive  ions,  the  more  complex  ions,  T(NaI)(HCONH2)m,  are  not 
detected  upon  dilution. 
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Mass-resolved  retarding  potential  measurements  are  shown  in  Fig.  4  for 
Na+(HCONH2)4  and  high  mass  (m/q  >  400  amu)  emitted  charged  particles  obtained  by 
operating  the  quadmpole  in  an  rf-only  high-pass  mode.  The  capillary  potential  was 
biased  at  550  V  while  the  retarding  grid  analyzer  was  at  ground  (see  Fig. 2  for  a 
schematic  of  the  experimental  setup),  From  these  measurements  we  conclude  that 
charged  droplets  are  formed  at  a  potential  -200  V  below  the  capillary  potential,  while  the 
solvated  ions  are  formed  with  a  broader  energy  distribution  centered  slightly  below  the 
droplet  appearance  energy,  with  a  minor  contribution  at  higher  energies.  The  distributions 
look  very  similar  for  ions  with  other  solvation  numbers.  Meanwhile,  the  energy 
distributions  do  not  change  noticeably  for  negative  ions  and  the  dilute  solution. 

Although  no  simulations  were  carried  out  for  KI  solutions,  we  found  it  useful  for 
comparative  purposes  to  also  investigate  mass  spectrometically  Kl/formamide  solutions. 
A  solution  at  the  same  weight  concentration  as  the  concentrated  Nal  solution  froze  or 
crystallized  at  the  capillary  orifice  and  no  measurements  could  be  carried  out.  We 
succeeded  in  creating  stable  cone-jets  with  1/4  volumetric  diluted  and  1/8  volumetric 
solutions  at  extraction  voltages  of  2100  V  and  2300  V,  respectively.  The  same 
corresponding  solvated  ions  were  observed,  and  the  abundances  are  listed  in  Tables  4  and 
5.  As  for  the  Na+  ions,  the  simple  solvated  K+  ions  have  distributions  that  are  not  greatly 
affected  by  dilution,  with  peak  at  n  =  5,  and  significantly  higher  abundances  at  larger 
solvation  numbers  than  observed  for  Na\  Complex  ions,  K+(KI)(HCONH2)m,  are 
generally  more  abundant  than  in  the  Nal  solutions,  and  are  observed  for  both  dilute 
solutions.  The  solvated  anion  distributions  show  the  same  trend  of  peaking  at  higher 
solvation  numbers  upon  dilution  as  observed  for  the  Nal  solutions.  Contrary  to  the  dilute 
Nal  solution,  complex  anions  are  also  observed  for  the  KI  solutions. 

4,2  Solvated  ion  thermochemistry  and  ion  survival  fractions 

The  binding  energies  of  the  calculated  solvated  ions,  Na+(NH2CHO)„,  n  =  3-6,  and 
['(NH2CHO),,,  n  =  2-5,  determined  from  the  quantum  chemical  calculations,  are  shown  in 
Table  6.  Also  listed  are  the  zeropoint  energies,  E:pe ,  and  the  internal  energies,  Ewl,  above 
the  zeropoint  energies  at  298.15  K  (difference  between  thermal  energy  and  zeropoint 
energy).  We  find  that  the  internal  energy  at  298.15  K  exceeds  the  binding  energy  above 


20 


n  =  4  for  the  cationic  solvated  ions,  and  n  —  3  for  the  1*  solvated  ions,  respectively.  Table 
6  also  lists  the  corresponding  classical  values  determined  for  the  clusters  bound  by  the 
interatomic  potentials  used  in  the  MD  simulation.  In  the  case  of  the  Na+  clusters,  the  MD 
binding  energies  are  significantly  lower  than  the  quantum  chemical  values  at  low  n  with 
improved  agreement  as  n  increases.  The  agreement  is  generally  good  between  the  MD 
binding  energies  and  the  quantum  chemical  values  for  the  T  clusters.  The  MD  internal 
energies  correspond  to  the  total  cluster  kinetic  and  potential  energy  minus  the  minimum 
potential  energy,  with  the  exclusion  of  the  kinetic  energy  associated  with  the  center-of- 
mass  velocity.  The  energies  are  obtained  from  averages  over  random  samples  of  clusters 
taken  from  separate  MD  simulations  of  such  cluster  ensembles  thermalized  initially  at 
310  K,  As  is  evident  that  the  onset  of  thermal  instability  is  for  clusters  with  the  number 
of  formamide  ligands  above  n  =  5  for  Na+  and  n  =  4  for  T  ions.  The  satisfactory 
agreement  with  the  quantum  chemical  predictions  supports  the  conclusions  that  we  drew 
from  comparison  between  the  experimental  and  the  MD  results  about  the  properties  of 
the  field-evaporated  ions. 

Table  7  lists  the  thermal  ion  survival  fractions  calculated  for  selected 
Na+(NHjCHO)„  ions  evaporated  with  characteristic  temperatures  of  300,  350  and  400  K. 
The  survival  fractions  are  calculated  for  infinite  observation  times,  mass  spectrometer 
flight  times,  and  the  MD  simulation  flight  time  of -15  ps,  and  are  based  on  the  RKK.M- 
type  unimolecular  dissociation  rate  calculations  described  in  Section  3.2.  The 
corresponding  survival  fractions  for  selected  TfNHaCHO),,  ions  are  given  in  Table  8. 

4.3  Molecular  dynamics  simulations 

Dielectric  and  charged  droplets  in  uniform  electric  fields  are  well-studied  problems 
involving  electrified  fluids.’1  25, 52-5 5  The  salient  features  of  the  physics  associated  with 
electrosprays  operating  in  the  cone-jet  mode  may  be  seen  already  in  the  simpler  problem 
of  a  conducting  liquid  droplet  elongated  in  a  uniform  external  electric  field,  including  the 
formation  of  electrified  jets,  as  well  as  the  emission  and  acceleration  of  charged  species 
at  the  ends  of  the  elongated  parent  drops.  Indeed,  Taylor2'  studied  the  instabilities  and  jet 
formation  of  fluids  emerging  from  an  electrospray  device  as  a  way  to  investigate  in  a 
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more  controlled  manner  the  physics  related  to  the  elongation  and  jet  formation  of  droplets 
in  a  uniform  electric  field. 

In  light  of  the  above,  atomistic  molecular  dynamics  (MD)  studies  were  initiated  of 
field-induced  fluid  instabilities  relevant  to  the  operation  of  electrospray  devices  by 
studying  first  the  closely  related  phenomena  present  in  the  behavior  of  dielectric  droplets, 
with  and  without  dissolved  salts.  As  discussed  below,  the  results  of  the  MD  simulations 
complement  and  agree  with  the  experimental  electrospray  resuits  to  a  high  degree.  First, 
we  present  results  from  the  study  of  pure  formamide  droplets,  of  the  same  size  as  the 
formamide/Nal  solution  droplets  that  are  our  primary  focus. 


4.3.1  Simulations  of  pure  formamide 

Since  the  potential  parameters  employed  in  our  simulations  were  not  specifically 
developed  to  reproduce  certain  physical  properties  of  the  simulated  formamide  (such  as 
the  melting  point  and  dielectric  constant)  it  was  necessary  to  perform  preparatory 
simulations  to  ascertain  the  physical  behavior  of  the  simulated  liquid  formamide. 
Constant  pressure  simulations  and  liquid-solid  coexistence  studies8’  indicate  that  our 
description  of  formamide  yields  a  liquid  density  that  is  to  within  a  few  percent  of  the 
experimental  value,  and  a  melting  point  that  is  approximately  10  K  above  the 
experimental  value  of  275  K.  For  an  improved  comparison  with  experiment,  and  to  avoid 
freezing  of  the  field-evaporating  salt  solution  droplets,  the  simulations  are.  therefore, 
performed  at  310  K  (that  is,  10  K  above  the  experimental  conditions). 

The  droplets  used  in  our  simulations  were  created  by  ‘carving  out’  a  spherical 
region  from  an  equilibrated  bulk  liquid  simulation  that  was  performed  at  constant 
temperature  and  zero  pressure,  with  a  subsequent  equilibration  prior  to  the  application  of 
an  external  field.  The  diameter  of  the  spherical  droplet  chosen  for  the  MD  simulations  is 
10  nm  and  it  contains  7150  formamide  molecules  (each  molecule  comprised  of  six 
atomic  partial  charges).  This  droplet  is  large  enough  to  allow  meaningful  comparison 
with  continuum  electrohydrodynamic  theories,  and  in  addition  it’s  curvature  is  similar  to 
the  estimated  local  curvature  of  the  experimentally  probed  liquid  jet  at  the  place  where 
the  electric  is  maximal  (see  Eqs.  6  and  7).  From  a  practical  point  of  view,  the  chosen 
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droplet  size  is  sufficiently  small  to  allow  a  wide  range  of  computer  experiments  to  be 
performed  in  a  reasonable  amount  of  time. 

The  droplet  was  equilibrated  in  a  range  of  increasing  electric  field  strengths  (from 
0  to  3  V/nm).  The  effect  of  the  externally  imposed  electric  field  can  be  seen  in  Fig.  5 
where  we  display  snapshots  of  the  droplet  for  three  values  of  the  applied  field,  illustrating 
instability  at  a  critical  electric  field  strength.  At  field  strengths  below  -0.5  V/nm,  the 
droplet  undergoes  mild  elongation  with  increasing  field  strength.  At  a  field  strength  of 
0.625  V/nm  the  droplet  becomes  unstable  and  undergoes  a  pronounced  elongation 
resulting  in  a  change  in  aspect  ratio  (r|=  c/a,  where  c  and  a  are  the  major  and  minor  axis 
of  the  ellipsoid)  from  q  -1  to  q  -13.  The  elongation  is  also  reflected  in  changes  in  the 
average  average  molecular  dipole  moment  projected  onto  the  z  axis.  There  have  been  a 
number  of  studies  addressing  the  critical  external  field  at  which  a  dielectric  droplet 
displays  this  instability.  In  early  work,  Taylor'5  theoretically  predicted  that  a  droplet  of 
initial  radius  R  (in  zero  external  field)  becomes  unstable  at  an  external  field  strength, 

E enticah  satisfying  (MKS  units): 

1 .625(r/2TO,rf)l,!,  (ID 

Applying  y  =  0.046  N/m  for  the  surface  tension  of  the  present  droplet  (computed  in 
separate  simulations)81  and  d=  9.68  nm  for  the  computed  equimolar  diameter,  results  in 
a  predicted  critical  applied  electric  field  of  E  =  0.48  V/nm  which  is  only  slightly  lower 
than  the  interval  of  0.5-0.625  V/nm  in  which  the  instability  is  observed  to  occur  in  the 
MD  simulation.  We  note  that  if  the  experimental  surface  tension82  y  =  0.0575  N/m  is 
used,  the  predicted  critical  field  strength  is  0.53  V/nm.  We  conclude  that  a  continuum 
formalism  provides  an  adequate  description  of  our  results  obtained  for  nanometer  scale 
droplets.  Additionally,  it  was  discovered81  that  the  external  electric  field  can  induce 
crystallization  of  the  pure  formamide  droplet. 

4.3.2  Simulations  of  10  nm  Nal-formamide  droplets 

Having  illustrated  that  the  thermodynamics  and  phase  behavior  of  a  dielectric  droplet  can 
be  changed  and  manipulated  by  external  electric  fields,  we  turn  next  to  investigations  of 
formamide  droplets  with  a  dissolved  Nal  salt  at  concentrations  that  are  comparable  with 
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those  used  in  the  experiments  described  in  section  4.! .  Simulations  were  performed  on 
10  nm  Nal-formamide  solutions  with  Nal/formamide  ratios  of  1/8  and  1/16,  respectively 
in  fields  of  0.5,  0.625,  0.75,  1.0  and  1.5  V/nm.  The  concentrated  droplets  contain  7150 
formamide  molecules  and  900  dissolved  Nal  units,  corresponding  to  a  29%  wt%  solution 
that  is  nearly  identical  to  the  concentrated  solution  used  in  the  experiments.  The  salt 
solution  droplets  were  formed  in  the  same  manner  as  the  pure  formamide  droplets 
described  earlier. 

Snapshots  of  the  1/16  salt  solution  droplet  subjected  to  an  electric  field  of  0.625 
V/nm  are  shown  in  Fig.  6.  A  stage  of  initial  elongation  along  the  field  direction  is  shown 
in  Fig  6b,  400  ps  after  the  field  has  been  turned  on.  Fig.  6c,  taken  300  ps  later,  exhibits 
further  elongation  of  the  droplet,  and  the  onset  of  charge  emission  in  the  form  of  a  Na* 
solvated  ion.  In  the  snapshot  shown  in  Fig.  6d,  which  was  recorded  1 . 1  nsec  after  the 
initial  application  of  the  field,  formation  of  jets  emanating  from  both  ends  of  the 
elongated  droplet  are  apparent.  While  a  large  solvated  ion  (with  about  20  formamide 
ligands)  with  a  single  Na+  ionic  at  the  center  is  seen  to  be  on  the  verge  of  being  emitted  at 
the  positive  end  of  the  parent  droplet  (on  the  right  in  Fig.  6e),  a  large  complex  negative 
ion.  r(NaI)(NIl2CHO)„,  is  seen  to  be  detaching  at  the  other  (negative)  end  of  the  droplet 
(left  in  Fig.  6e).  About  25  ps  after  the  instant  when  the  configuration  shown  in  Fig,  6d 
was  recorded  ,  a  section  of  the  material  streaming  from  the  positively  charged  side  of  the 
droplet  thins  down  to  a  narrow  hydrogen  bonded  filament  of  formamide  molecules 
resulting  in  the  detachment  of  a  cluster,  or  nanodroplet,  of  charge  +3  (see  Fig.  6e) 
consisting  of  17  Na"  ions  ,14  T  ions  and  -200  formamide  molecules.  This  corresponds 
to  a  mass-to-charge  ratio  of  about  3,700  Daltons. 

The  charged  clusters  that  detach  from  the  two  ends  of  the  droplet  are  accelerated 
by  the  applied  electric  field  and  their  physical  characteristics  (charge,  mass,  size, 
energies),  are  recorded  when  they  cross  two  planes  (which  we  refer  to  as  "detection 
planes)  located  at  ±  200  nm  from  the  droplet  center  along  the  field  direction  (the  z-axis). 
Properties  of  the  recorded  positively  and  negatively  charged  clusters  are  summarized  in 
Fig.  7  as  a  function  of  the  clusters’  radii,  estimated  by  using  or  each  cluster  the 
equimolar  radius  defined  by  the  cluster  mass,  m,  the  computed  density  of  the 
Nal/formamide  solution,  p-1.32  g/cm3,  and  the  relation  m  =  4/37tR"1p-  Figure  7a  displays 
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the  estimated  surface  electric  field,  E,  of  each  cluster  using  Coulombs  law  (Eq.  2).  A 
change  in  the  dependence  of  E  on  the  cluster  radius  R  is  observed  at  R  «  1  nm,  with 
cluster  having  radii  R  <  1  nm  corresponding  to  single  ions  solvated  by  formamide 
molecules,  while  clusters  with  R  >  1  nm  are  characterized  as  multi-ion  droplets,  In  the 
latter  dispersion  forces  may  be  expected  to  make  a  relatively,  together  with  the 
electrostatic  interactions.  The  change  in  behavior  occurring  for  clusters  with  R  >  1  nm  is 
reflected  also  in  the  q/m  and  q/qRayieigh  ratios  plotted  in  Fig.  7b  and  7c,  respectively.  The 
vertical  dashed  line  in  these  figures  corresponds  to  the  experimental  estimate  (see  Table 
1)  for  the  radius  curvature  of  the  electrospray,  Rmax,  at  the  point  of  maximum  electric 
field,  Emax  (see  Eqs.  6  and  7).  Essentially  all  the  clusters  recorded  in  the  MD  simulation 
have  radii  smaller  than  Rm&.  The  few  clusters  with  radii  larger  than  /fmax  appeared  at  the 
very  end  of  the  simulation  when  breakup  of  the  droplet  into  large  fragments  was 
observed.  The  experimental  estimate  for  Em3X  =  1 .67  V/nm  (  determined  from  Eqs.  6  and 
7  as  the  a  average  between  the  positive  and  negative  polarities)  is  indicated  by  the  thin 
solid  line  in  Fig  7a.  The  estimated  surface  field  of  the  simulated  clusters  remains  very 
close  to  /imax  for  cluster  radii  in  the  droplet  regime. 

Figure  7b  shows  the  charge-to-mass  ratios  for  the  clusters  emitted  by  the 
simulated  droplet.  The  mean  value  of  the  experimental  q/m  (Table  1 )  is  shown  as  the 
horizontal  dashed  line.  The  MD  values  quickly  approach  the  mean  experimental  value  as 
the  cluster  sizes  enter  the  droplet  regime.  If  one  uses  Emn  as  a  rough  limiting  value  of  the 
surface  field  for  a  charged  droplet  in  Eq.  (2)  (Coulomb’s  law)  to  write  i/=47rR'£o£max,  an 
estimate  for  the  upper  bound  of  q/m  can  be  obtained  from; 

q/m  =  3e0Emm/pR. 

The  upper  bound  for  q/m  as  a  function  of  £mJW  and  R  is  is  shown  in  Fig.  7b  as  a  thin  solid 
line.  It  is  evident  that  the  data  for  charged  droplets  (R  >  1  nm  recorded  in  the  MD 
simulations,  follow  this  line. 

Figure  7c  shows  the  ratio  between  the  observed  (MD)  charge  and  the  predicted 
charge  for  Rayleigh  instability  and  fission,  computed  from  Eq.  (1),  where  the  surface 
tension  y  =  0.046  N/m  determined  from  the  MD  formamide  droplet  is  used  as  an 
approximation  for  the  salt  solution.  For  clusters  in  the  droplet  regime,  their  charge 
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exhibits  variations  about  0.5  of  the  Rayleigh  limit.  Using  Eq.  (1)  and  (2)  one  obtains  the 
upper  limit  for  the  ratio.  qlqRayieigh- 
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The  corresponding  values  of  Eq.  (13)  are  shown  as  a  solid  line  in  Fig.  7c.  This  simple 
estimate  captures  the  general  trend  and  values  observed  in  the  MD  simulation. 

Fig.  8  compares  profiles  of  ion  and  solvent  molecule  densities  along  the  semi- 
major  axis  of  the  droplet  spheroids  for  the  cases  shown  in  Figs.  6  (a-c).  The  figure  also 
provides  profiles  of  the  local  electric  field  and  the  charge,  the  latter  given  by  the 
difference  between  positive  and  negative  ion  densities.  The  plotted  electric  field  does  not 
include  the  external  field  (0.625  V/nm).  The  densities  are  calculated  by  defining  prolate 
spheroid  shells  that  match  the  instantaneous  droplet  shape  and  end  curvatures,  and  by 
counting  the  respective  species  between  neighboring  shells  within  a  1  nm  core  radius 
with  respect  to  the  z-axis.  While  the  radial  diameter  of  the  narrowing  droplet  ends  are  less 
than  that  of  the  sampled  core  diameter  (2  nm),  the  statistics  of  the  density  derive  only 
from  material  inside  the  prolate  spheroid  shells,  thereby  yielding  an  accuracy  superior  to 
that  obtained  via  simple  binning  of  the  z-axis.  This  process  was  not  implemented  for  the 
case  of  the  droplets  in  Figures  6d  and  6e  due  to  the  narrow,  rapidly  varying,  curvature. 
The  data  are  obtained  from  300  ps  averages  in  the  vicinity  of  the  snapshots  shown  in  Figs 
6  (a-c). 

The  electric  field  strength  along  the  z-axis  is  computed  by  binning,  (during  the 
simulation)  the  number  of  each  species  (Na,  I,  and  6  formamide  atoms  each  having  their 
own  partial  charge)  with  the  z-  and  r-  axis  having  bin  widths  of  0.2  nm.  The  binned  data 
is  averaged  over  a  specified  time  in  which  the  droplet  does  not  change  its  shape  or  size 
significantly.  Based  on  the  charge  of  each  species,  the  binned  number  of  atoms  provide 
the  average  amount  of  charge  coming  from  each  species  in  each  of  the  radial  'rings' 
associated  with  each  z-bin.  The  potential  V(Z)  at  a  point,  Z,  on  the  z-axis  is  then  a  sum 
of  the  contributions  from  the  total  charge  q,  in  each  of  the  binning  rings  of  radius  r\  for 
all  the  z-bins.  (except  for  the  central  radial  bin  in  the  z-bin  for  which  the  field  is 
determined): 
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The  electric  field  at  point  Z  without  inclusion  of  the  external  field  is  then  given  by 

E(Z)=-dV(Z)/dZ. 

The  field-free  droplet  exhibits  a  depletion  of  ions  at  the  interface,  with  a  slightly 
higher  ion  density  in  the  second  or  third  interfacial  layer  (-6  A  from  the  interface)  in 
comparison  with  the  bulk.  Differences  between  Na"  and  1'  densities  are  insignificant. 
This  is  consistent  with  the  traditional  view35  pertaining  to  the  behavior  of  the  interfacial 
region  of  salt  solutions,  and  it  contrasts  with  recent  MD  and  experimental  findings  in 
aqueous  Nal  solutions,  where  the  highest  T  densities  were  found  in  the  outer  interfacial 
layer,3*'311’47’49;  we  defer  further  comments  on  this  issue  to  the  discussion  section. 

Figure  8b  shows  the  profiles  for  the  droplet  shown  in  Fig.  6b.  The  interfacial 
charge  density  has  clearly  increased.  In  case  of  the  positive  end  (right  in  Fig.  6),  there  is 
a  slight  preference  for  positive  ions  with  ion  densities  of  both  polarities  peaked  at  the 
interfacial  layer.  At  the  negative  droplet  end,  there  is  a  clear  preference  for  negative  ions 
in  the  interfacial  layer.  The  electric  field  strength  is  highest  beyond  the  droplet  ends, 
with  total  fields  (including  the  external  field)  exceeding  1.5  V/nm.  Interestingly,  the 
maxima  are  quite  far  from  the  ion  density  maxima,  suggesting  that  the  aligned  dipoles  of 
the  solvent  molecules  contribute  considerably  to  the  field  right  outside  the  droplet. 
Indeed,  when  the  field  produced  by  the  droplets  is  decomposed  into  contributions  from 
the  ions  and  the  formamide  molecules,  we  find  that  they  each  make  a  significant 
contribution.  It  is  also  of  interest  to  note  that  the  axial  electric  field  of  the  pure 
formamide  droplet  (placed  in  an  external  field  of  0.5  V/nm,  see  Fig  5)  is  essentially  the 
same  as  that  shown  in  Fig.  8b  for  the  ionic  solution  droplet.  In  the  pure  formamide 
droplet  the  total  contribution  to  the  internal  field  comes  from  the  aligned  diploes  of  the 
formamide  molecules.  From  the  fact  that  the  internal  fields  are  of  similar  values  for  the 
pure  and  the  salt  solutions  drops,  we  conclude  that  the  larger  contribution  of  the 
molecular  dipoles  field  in  the  pure  drops  reflects  a  higher  degree  of  alignment  in  these 
drops  compared  to  the  salt  solution  drops  where  the  dissolved  ions  frustrate  the  molecular 
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alignment  in  the  field  direction.  Direct  analysis  of  the  dipole  alignments  in  pure 
formamide  and  salt  solution  drops  verifies  the  above  conclusion.  We  also  note  that  the 
external  field  strengths  of  the  order  of  V/nm  originating  from  droplet  polarization  are 
established  very  shortly  subsequent  to  application  of  the  the  external  and  before  the 
occurrence  of  any  significant  field-induced  elongation. 

The  field  produced  by  the  droplet  vanishes  asymptotically  far  from  the  droplet  so 
that  the  total  summed  field  converges  to  the  external  constant  field  of  0.625  V/nm.  The 
negative  field  inside  the  droplet  is  of  the  same  nature  as  the  field  within  a  dielectric 
slab  or  sphere  subjected  to  an  external  field.  Charges  of  opposite  polarity  are  induced  on 
opposing  sides  of  the  dielectric  and  this  creates  an  internal  field  that  opposes  the  external 
field. 

In  the  present  case,  the  average  internally  induced  field  is  —0.58  V/nm, 
corresponding  to  a  small  total  internal  field,  E  =  0.625  -  0.58  -  0.04  V/nm.  Such  smalt 
internal  fields  are  typical  of  dielectric  materials  responding  to  external  electric  fields. 
The  high  sensitivity  of  computed  dielectric  properties  to  small  variations  in  this  internal 
field  necessitates  the  accumulation  of  sufficient  data  that  would  permit  meaningful 
statistical  analysis.  The  highly  oscillatory  nature  of  the  field  within  the  drop  seen  in  our 
simualtions,  is  due  to  limited  statistics  and  the  fact  that  the  computed  field  reflects  a  non¬ 
uniformity  of  the  ionic  distribution  that  does  not  self-average  over  the  averaging  interval 
of  -300  ps  used  in  these  simulations. 

4.3.3  Charge  emission  statistics 

Statistics  on  charge  emission  is  gathered  by  recording  the  passage  of  all  ejected  clusters 
through  the  aforementioned  “detection  planes”  placed  at  ±200  nin  along  the  z-axis.  The 
atomic  position  and  momentum  data  of  each  cluster  is  also  recorded  so  that  the  evolution 
of  an  individual  cluster  can  be  examined  in  a  separate  continued  simulation  of  the 
isolated  species. 

The  surface  tension  of  the  Nal -formamide  droplets  is  comparable  to  that  of  the  pure 
formamide  droplet  and  the  predicted  critical  field  for  instability25  discussed  earlier  {  see 
also  Fig.  5)  is  similar.  However,  in  the  salt-solution  case,  the  droplet  does  not  simply 
elongate  to  another  stable  configuration  at  the  critical  field.  Instead  the  imbalance 
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between  electrical  stresses  and  capillary  forces  is  relieved  through  a  more  favorable 
route,  namely  one  involving  formation  of  a  jet  of  charged  clusters  (at  each  of  the 
elongated  nano-drop).  At  the  smaller  investigated  fields  (<  1  V/nm),  the  droplets 
elongate  more  slowly  and  emit  cluster  ions  less  frequently  with  the  result  that  capillary, 
Rayleigh-type,  instabilities  lead  to  break-up  of  the  droplet  into  large  charged  segments 
that  are  accelerated  out  of  the  computational  cell  before  the  number  of  cluster  ions  that 
are  emitted  are  adequate  for  meaningful  statistics.  Consequently,  when  investigating 
cluster-ion  emission  statistics  we  focus  on  simulations  of  drops  at  the  higher  field 
strengths.  The  improved  statistics  associated  with  the  higher  charge  fluxes  provide  a 
better  comparison  with  the  mass  spectrometric  measurements. 

Under  the  influence  of  an  electric  field  of  1.5  V/nm,  the  droplet  begins  to  deform 
about  10  to  20  ps  prior  to  charge  emission,  along  the  axial  field  direction,  of  positively 
and  negatively  charged  solvated  ions  from  opposite  ends  of  the  drop.  Over  the  course  of 
about  1  ns,  the  parent  drop  continues  to  elongate  to  at  most  an  aspect  ratio  q  =  8,  and  the 
onset  of  charge  emission  is  observed.  The  continuous  stream  of  charged  particles 
prevents  elongation  to  higher  aspect  ratios  such  as  the  one  observed  for  the  pure 
formamide  droplet  (see  Fig.  5).  The  simulations  demonstrate  that  charge  emission 
happens  from  multiple  sites  on  each  end  of  the  drop.  This  is  consistent  with  a  stressed, 
multi-cone  electrospray  regime  discussed  earlier.  We  remark  here  that  because  of  the 
early  emission  of  charged  particles  in  simulations  employing  a  high  electric  field,  ( 1 .5 
V/nm)  we  focused  on  lower  applied  fields  which  enable  comparisons  with  the 
experiments,  particularly  of  issues  pertaining  to  the  dynamics  of  the  particle  emission  and 
droplet  formation. 

The  charge  distribution  of  cluster  ions  originating  from  formamide  droplets 
containing  29  %  (by  weight)  of  dissolved  Nal  is  shown  in  Fig.  9.  Results  are  presented 
here  for  two  separate  simulations  having  externally  applied  electric  field  strengths  of  1 .0 
and  1.5  V/nm.  The  abundances  of  charge  of  both  these  studies  are  similar  showing  that 
most  of  the  observed  clusters  have  charge  ±1,  however,  multiply  charged  species  are  also 
observed. 

A  more  detailed  information  on  the  ionic  composition  of  the  most  frequently 
observed  clusters  of  charge  ±1  is  shown  in  Fig.  10.  In  case  of  singly  positively  charged 
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clusters,  the  fraction  of  clusters  with  ionic  composition  Na+,  Na2+I",  etc,  and  for 
negatively  charged  clusters,  the  fraction  corresponding  to  I",  Na+Ij ",  etc.  are  shown 
without  regard  to  the  number  of  accompanying  adsorbed  formamide  molecules.  While 
multi-ion  clusters  contribute,  the  primary  contribution  is  from  single-ion  clusters  (Na+ 
and  r  with  varying  numbers  of  adsorbed  formamide  molecules).  Furthermore,  the  multi¬ 
ion  cluster  fraction  increases  with  field.  This  is  consistent  with  the  present  experimental 
single  cone-jet  observations  at  lower  representative  field  strengths,  where  the  multi-ion 
fraction  due  to  Na2+I~  and  Na+l2~  solvated  ions,  is  very  small  (see  Fig.  3  and  Tables.  1-4). 

In  Fig.  1 1  we  compare  the  solvation  number  distribution  of  the  singly-charged 
positive  and  negative  cluster  ions  compiled  from  the  concentrated,  1,5  V/nm  MD 
simulation  to  those  observed  in  the  concentrated  electrospray  experiment.  Also  shown 
are  MD  results  adjusted  for  a  -30  ps  flight  time  using  the  calculated  survival  fractions 
from  Tables  7  and  8  for  300K.  For  Na‘  and  V  solvation  numbers  exceeding  6  and  5, 
respectively,  a  zero  survival  fraction  is  assumed,  while  for  n  <  4  and  3,  respectively,  a 
!  00  percent  survival  is  implied.  The  Na+  clusters  have  predominantly  4-5  formamides, 
and  the  adjusted  cluster  distribution  compares  well  to  the  experimental  abundances.  The 
I’  clusters  exhibit  a  much  broader  distribution  with  a  bias  toward  lower  solvation 
numbers.  The  adjusted  MD  statistics  are  again  consistent  with  the  experimental  results 
(Tables  2  and  3)  for  a  single  cone-jet  electrospray  of  the  same  concentrated  solution. 

5.  Discussion 

The  close  correspondence  that  we  found  between  the  data  measured  in  the  vacuum 
electrospray  experiments  and  the  results  recorded  in  our  theoretical  MD  simulations 
performed  on  liquid  drops  (made  of  pure  formamide  or  Nal/formamide  solutions), 
demonstrate  that  such  coordinated  , complementary,  and  supplementary  studies  are  most 
valuable  sources  of  information  and  deep  insights  pertaining  to  the  physics  and 
operational  control  of  electrospray,  nanojets,  and  charged  droplets  under  the  influence  of 
high  electric  fields.  Having  presented  in  the  previous  sections  our  laboratory 
measurements  and  theoretical  simulations  we  discuss  now  certain  issues  that  emerged 
through  the  analysis  of  the  results  of  our  investigations. 
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The  characteristic  parameters  of  our  experiments  (see  Table  I)  demonstrate  that 
the  cone-jet  generated  in  these  experiments  operates  in  the  nanoscale  jet  regime  with 
sufficiently  high  fields  generated  normal  to  the  liquid  surface  to  induce  field  evaporation 
of  ions.  The  observed  average  mass-to-charge  ratios  and  the  significant  cluster-ion 
currents  are  signatures  of  a  mixed  droplet/ion  emission  regime.  The  presence  of  field- 
induced  ion  evaporation  is  consistent  with  the  estimated  values  for  the  maximum  surface- 
normal  fields  Emax  (between  1.4  and  1.7  V/nm)  and  for  the  critical  radii,  Rmax  (between  2 
and  3  nm)  derived  from  the  axial  emission  currents  and  mass  flows,  thus  confirms 
predictions  obtained  from  the  theory  of  cone-jet  emission.  I5'"  The  derived  jet  diameter 
of  about  5  nm  is  comparable  to  that  of  the  simulated  droplet,  supporting  the  comparison 
between  the  present  electrospray  experiment  and  the  MD  droplet  simulations. 

Decreasing  the  salt  concentration  of  the  Nal/formamide  solution  leads  to  a 
reduction  of  the  maximum  surface-normal  field  strength  and  an  increase  of  the  associated 
radius  of  curvature.  This  is  consistent  with  Eqs.  (6)  and  (7).  The  expected  increase  in 
m/q,  however,  is  found  to  be  significant  only  for  negatively  charged  emitted  particles. 
This  somewhat  surprising  result  can  be  explained  by  the  substantial  decrease  in  charged 
particle  beam  divergence  upon  dilution,  as  reflected  by  the  ratio  of  the  Faraday  Cup 
current  with  respect  to  the  total  emitted  current.  As  angular  distribution  measurements  on 
ionic  liquid  cone-jets  operated  in  a  nano  regime  have  shown,  "  the  m/q  distribution  can 
be  non-uniform  with  respect  to  the  emission  polar  angle.  Thus,  the  derived  average  m/q 
value  is  only  representative  of  charge  emission  captured  at  small  angles  with  respect  to 
the  capillary  axis.  The  non-uni formity  reported  by  Chiu  et  al.  for  the  ionic  liquid, 
[Emim][Im],  ’  is  consistent  with  predictions  that  the  surface  electric  field  is  highest  at 
the  cone-jet  neck,  where  field-evaporated  ions  are  driven  away  from  the  cone-jet  axis."' 
31  The  respective  field-evaporated  ions  are,  therefore,  preferably  emitted  at  large  angles 
while  droplets,  emitted  from  the  jet  tip,  have  narrow  angular  distributions  centered  on 
axis.  This  interpretation  was  further  confirmed  for  the  ionic  liquid  by  the  measured 
energy  distributions  demonstrating  that  the  ions  peak  at  energies  nearly  equal  to  the 
emitter  potential,  while  the  droplets  appear  at  lower  potentials.  Consequently,  if  similar 
field  evaporation  dynamics  apply  in  the  present  concentrated  Nal/formamide  solution,  a 
significant  fraction  of  the  ion  current  is  not  detected  by  the  Faraday  Cup,  and  the  average 
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m/q  of  the  entire  jet  is  considerably  higher  for  cases  with  substantial  divergence  of  the 
ionic  component  in  comparison  with  the  droplet  component.  This  interpretation, 
however,  is  contradicted  by  the  fact  that  the  present  retardation  potential  measurements 
(Fig.  4)  show  that  only  a  small  fraction  of  ions  are  produced  at  energies  exceeding  the 
droplet  energy,  suggesting  that  most  ions  are  formed  in  a  region  of  the  neck  downstream 
from  where  ohmic  losses  have  occurred.  Similar  retarding  potential  results  have  been 
derived  by  Gamero-Castano  for  concentrated  Nal/formamide  at  low  flow  rates.24 
Consequently,  the  field  evaporation  dynamics  of  the  present  solution  may  not  be 
equivalent  to  that  observed  for  the  ionic  liquid  [Emim][Im]  This  does  not  discount  the 
possibility  that  ion  emission  is  more  divergent  than  the  droplet  spray. 

In  addition  to  field  evaporation  from  the  jet,  ions  can  be  emitted  from  droplets  as 
described  by  the  mechanism  of  Iribame  and  Thomson6  for  electrosprays  of  volatile 
solvents  in  an  ambient  atmosphere.  This  requires  that  the  surface  electric  field  strength, 
given  by  Eq,  (2),  exceeds  the  ion  evaporation  limit.  However,  as  pointed  out  by  Gamero- 
Castano  and  Fernandez  de  la  Mora,22  this  mechanism  can  only  lead  to  a  current  that  is  a 
small  fraction  of  the  total  droplet  current  because  of  the  fast  decrease  of  the  electric  field 
in  the  direction  normal  to  the  liquid  surface  following  the  emission  of  ions. 

Further  insights  were  obtained  through  comparative  mass  spectrometric  studies  of 
concentrated  and  dilute  Nal/formamide  solutions.  While  the  distribution  of  Na*  solvated 
ions  is  very  similar  at  the  two  concentrations,  apart  from  a  greater  number  of  trace  ions  at 
high  solvation  numbers,  the  1"  solvated  ion  distribution  broadens  considerably  and  shifts 
to  higher  solvation  numbers  upon  dilution.  This  behavior  is  consistent  with  recent 
theoretical  and  experimental  investigations  of  aqueous  salt  solutions,34'40'  43-49  where  a 
preference  at  the  interface  for  more  polarizable  anions  such  as  Y  was  explained  in  part 
with  the  broader  range  of  solvation  numbers  that  the  anion  can  accommodate  in 
comparison  with  less  polarizable  anions  and  cations,  which  prefer  more  rigid,  hydrogen- 
bonded  solvation  shells.  Thus,  as  the  concentration  decreases,  more  solvent  molecules 
are  available  to  bond  with  F,  while  Na4  maintains  a  tightly  bound  structured  solvation 
shell. 

The  Na+  solvation  number  (the  number  of  solvating  formamide  (FA)  molecules 
per  solvated  Na+,  Na+(FA)n)  drops  rapidly  for  solvated  ions  with  n  >  4,  for  which  room 
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temperature  thermal  instability  is  predicted.  However,  an  increased  number  of  ions  with 
high,  thermally  unstable,  solvation  numbers  is  observed  upon  dilution.  This  suggests  that 
a  significant  number  of  cluster  ions  are  emitted  with  large,  unstable,  solvation  numbers, 
and  subsequently  they  lose  (via  evaporation)  some  of  the  solvating  ligands,  resulting  in 
evaporative  cooling  of  the  product  cluster  ions.  This  effect  is  even  more  extreme  for  the 
f  cluster  ions  to  the  point  that  the  most  abundant  solvation  number  lies  beyond  the 
present  predictions  of  thermal  stability  at  a  liquid  temperature  of  298  K  (Table.  6). 

Similar  behavior  was  observed  in  our  experiments  on  KI  solutions,  where  the  K+ 
solvated  ions  also  have  maximum  abundances  near  the  predicted  limit  of  thermal 
stability.  The  main  difference  between  the  KI  and  Nal  solutions  are  the  survival,  in  the 
KI  case,  of  complex  ions  (of  both  polarities)  upon  dilution  for  both  ion  polarities.  The  KI 
experiments  were  mired  by  crystallization  of  the  solute  and,  therefore,  had  to  be 
conducted  at  significantly  lower  molar  concentrations.  The  higher  propensity  for  forming 
complex  ions,  X^KIXHCONHjJm,  despite  the  lower  concentration  in  comparison  with 
the  Nal  solution,  may  be  a  manifestation  of  the  lower  solubility  of  the  KI  solute. 

The  snapshots  of  the  drop  configurations  recorded  in  the  simulation  of  a 
Nal/Formamide  solution  droplet  with  a  (Nal):  16  (formamide)  ratio,  under  the  influence 
of  an  electric  field  of  0.625  V/nm,  reveal  a  number  of  features  that  are  in  direct 
correspondence  with  the  experimental  observations.  In  particular,  in  both  cases  we 
observe  mixed  droplet  formation  and  field  ion  evaporation  events,  Formation  of  nanojets 
is  clearly  observed  in  the  simulations  (See  Fig.  6d  where  a  jet  with  a  diameter  of  about 
2nm  is  shown  ).  However,  we  remark  that  the  jet  emerging  in  our  simulation  is  a  transient 
state  whose  further  development  prior  to  particle  emission  and  breakup  into  droplets  may 
be  limited  by  the  finite  size  of  the  simulated  drop.  In  the  experiment,  on  the  other  hand,  a 
continuous  (unlimited)  source  of  liquid  is  present,  thus  allowing  the  development  of 
steady  state  jet  (flow)  conditions.  Therefore,  in  order  to  allow  detailed  theoretical  study 
of  the  breakup  processes  it  is  desirable  to  slow  down  the  breakup  processes.  To  this  aim 
we  used  MD  drop  simulations  performed  with  lower  applied  electric  fields.  Interestingly, 
in  the  simulations  the  majority  of  cluster  ions  are  found  to  be  emitted  from  the  tip  of  the 
jet  rather  than  from  the  neck  region  next  to  the  Taylor  cone  (see  Fig.  1 ).  This  finding  is 
consistent  with  the  experimental  energy  distributions  of  solvated  ions  that  demonstrates 
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that  indeed  most  ions  are  emitted  closer  to  the  jet  region,  that  is,  beyond  the  cone-neck 
region  where  ohmic  losses  occur.  This  observation,  contrasts  continuum 
electrohydrodynamic  theory  predictions  of  maximum  surface  electric  field  strengths  at 
the  neck.22'2'1 

Because  of  the  small  dimensions  of  the  jets  formed  in  the  MD  simulations,  and 
deduced  from  analysis  of  the  experimental  measurements,  one  may  indeed  expect  that  a 
continuum  approach  would  fail  to  provide  an  adequate  description  of  the  evolution  and 
breakup  processes  of  the  complex  liquid  system  on  which  we  focus  here.  From 
inspection  of  the  drop  configurations  displayed  in  Fig.  6  we  observe  the  occurrence  of 
spatial  and  temporal  fluctuations  of  the  ion  density  that  are  expected  to  influence  in  a 
significant  manner  the  electric  field  distribution  and  the  consequent  stability  and 
evaporation/fragmetation  dynamics  of  the  system  under  study.  Indeed,  we  have  noted 
that  droplet  formation  events  (location  and  timing)  correlate  with  instances  of  ion 
bunching  in  the  jet  region.  These  observations  are  likely  to  be  of  particular  relevance  to 
ionic  solutions  where  bunching  of  the  solute  ions  (sodium  and  iodine  ions  in  our  case) 
occurs;  indeed,  in  the  case  of  an  ionic  liquid  (where  more  uniform  ion  distributions  are 
expected)  characteristics  that  indicate  emission  from  the  neck  region  have  been 
observed,  ’  even  though  the  jet  diameter  in  that  system  is  estimated  to  be  similar  to 
that  found  by  us  in  the  present  study. 

Additionally,  stochastic  thermal  effects  play  a  much  increased  role  for  fluid  jets  of 
nano-scale  dimensions,5 1  modifying  the  kinetics  of  fluid  pinch-off  from  the  predictions  of 
a  purely  deterministic  continuum  approach. 

An  important  prediction  obtained  from  the  MD  simulations  concerns  a  crossover 
as  a  function  of  the  size  of  the  emitted  species,  that  distinguishes  between  solvated 
single-ion  clusters  with  radii  R  <  1  nm,  and  multi-ion  droplets  that  emerge  for  larger  size 
fragments,  R  >  1  nm.  Such  crossover  is  indeed  found  experimentally,  and  it  may  be 
explained  through  electrostatic  considerations  (see  Fig.  7).  This  crossover  reflects 
stronger  bonding  in  the  smaller,  single-ion  solvated  clusters  due  to  direct  ion-dipole 
interactions  which  correlate  with  the  ability  to  sustain  larger  values  of  charge-to-mass 
ratios  compared  to  larger  size  clusters  (see  Fig.  7b). 
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Emission  of  clusters  that  include  ion  pairs  was  detected  in  both  the  experiments 
and  simulations.  While  only  species  containing  a  single  ion  pair  were  observed  in  the 
experiments  for  concentrated  Nal  formamide  solutions,  and  none  for  the  dilute  solution, 
higher  fractions  of  such  ions  with  species  containing  2  and  more  ion  pairs  were  found  in 
the  simulation  (Fig.  10).  Note  that  the  respective  radii  of  the  complex  ions  containing  ion 
pairs  exceed  1  nm,  and  exhibit  electrostatic  properties  of  droplets  according  to  the 
analysis  of  Fig.  7.  We  have  not  conducted  an  RRKM  analysis  and  thermo-chemical 
analysis  for  these  species  to  assess  whether  such  complex  ions  decay  faster  than  the 
experimental  time  scale.  Similarly,  no  multiply  charged  ions  were  detected  in  the 
experiments,  while  they  were  found  in  the  simulation  (Fig.  9).  I  ligher  field-strengths  do 
not  appear  to  be  the  source  of  these  findings  since  an  increase  of  the  field  in  the 
simulations  from  1  V/nm  to  1.5  V/nm  did  not  affect  the  observed  trends  (see  Figures  9 
and  10).  We  note  however,  it  can  be  said  that  earlier  vacuum  electrospray  studies 
conducted  with  glycerol  solutions  observed  a  multitude  of  solvated  ions  including 

7  17  1  "1 

metastable  ions  and  multiply  charged  species.  "  '  ‘ 

In  general  the  emitted  cluster  ions  recorded  in  the  MD  simulations  contain  more 
nascent  formamide  molecules  (that  is  they  are  characterized  by  larger  solvation  numbers) 
than  the  clusters  observed  in  the  experiments  where  typical  flight  times  (that  is,  the 
elapsed  time  between  the  instant  of  cluster  generation  and  the  moment  of  it’s  detection) 
are  of  the  order  of  30  psec.  The  recorded  occurrence  of  larger  solvated  ion  clusters  is 
likely  to  be  correlated  with  the  propensity  of  the  solvent  (formamide)  to  form  hydrogen 
bonds  (see  Figs.  6(c-e)).  Since  the  largest  solvated  ions  observed  in  the  MD  simulation 
are  thermally  unstable,  evaporative  loss  of  formamide  molecules  over  longer  flight  times 
can  effectively  lower  the  average  solvation  numbers  to  values  close  to  those  found 
experimentally  (see  Fig.  1 1 ).  The  observation  of  larger  solvated  ions  in  the  simulations 
suggests  that  such  larger  clusters  ions  are  indeed  generated  in  the  laboratory  experiment, 
and  that  subsequently  these  large,  thermally  (room  temperature)  unstable,  clusters  reduce 
their  size  via  a  cascades  of  evaporative  cooling  processes  that  occur  along  the  clusters’ 
trajectories  to  the  detector.  Indeed,  separate  simulations  on  large  thermally  unstable 
cluster  ions  show  that  the  loss  of  a  formamide  molecule  leads  to  an  approximate 
reduction  in  cluster-ion  temperature  of  50  K.  As  evident  from  the  data  shown  in  Table  7, 
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a  50  K  change  in  temperature  has  a  rather  significant  effect  on  the  stability  of  the  cluster 
ions.  We  remark  here  that 

the  MD  simulation  was  conducted  at  a  relatively  high  external  field,  corresponding  to 
stressed  conditions,  which  could  explain  the  increase  of  the  nascent  cluster-ion 
temperature  from  3 1 0  K  to  about  400  K  from  the  beginning  to  the  end  of  the  simulation. 
While  the  comparison  with  the  present  experimental  solvated  Na+  ion  distributions, 
predicted  thermal  stabilities,  and  RRKM  decay  rates  suggest  a  surface  temperature  below 
350  K,  a  200V  drop  between  the  capillary  potential  and  the  observed  ion  appearance 
potential  indicates  the  occurence  of  ohmic  heating  in  the  neck  region  of  the  cone-jet. 
Thus,  a  slightly  elevated  temperature  at  the  point  of  ion  field  evaporation  could  be 
expected, 

A  rather  convincing  illustration  that  the  emitted  charged  droplets  recorded  in  the 
simulations  are  consistent  with  the  values  deduced  from  the  experimental  data  for  the 
characteristic  maximal  electric  field,  Emax,  and  the  corresponding  local  curvature 
(determined  by  the  local  radius,  Rmax,  of  the  flowing  fluid  structure),  can  be  found  in 
Fig.7.  Except  for  droplets  emitted  at  the  very  end  of  the  simulation,  where  the  main 
droplet  disintegrates,  the  radii  of  all  the  emitted  droplets  are  found  to  be  smaller  than  the 
polarity  averaged  experimentally  determined  value  of  Rmax.  We  also  find  that,  coincident 
with  the  above,  the  surface  fields  of  the  droplets  are  close  to  the  experimental,  polarity 
averaged,  value  of  EmtJX,  which  is  found  to  be  close  in  value  to  the  interfacial  field 
determined  from  the  simulations  (see  Figs.  8b  and  8c).  Additionally,  the  charge-to-mass 
ratio  associated  with  the  experimental  maximal  value  of  the  electric  field,  Emax  (Eq.  1 2), 
is  found  to  agree  well  with  the  values  derived  for  the  emitted  charged  droplets.  Finally, 
the  observation  that  the  emitted  droplets  are  charged  to  within  0.4  to  0.5  of  the  Rayleigh 
limit  is  found  to  be  consistent  with  Eq.  (13)  when  using  the  experimental  value  of  Em]I. 

It  is  pertinent  to  remark  here  on  the  small  asymmetry  between  the  calculated 
negative  and  positive  ion  interfacial  density  profiles  (see  Fig.  8),  particularly  in  light  of 
recent  interest  in  the  interfacial  properties  of  aqueous  solutions, J4~49  Upon  application  of 
an  electric  field  the  ion  density  at  the  interface  increases  (see  Fig.  8b).  While  at  field-free 
conditions  no  discemable  difference  is  found  between  the  negative  and  positive  ion 
interfacial  densities,  a  simulated  drop  placed  in  a  uniform  electric  field  exhibits  a  small 
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preference  for  negative  ions  near  the  negative  polarity  interface,  while  positive  ions  are 
only  marginally  preferred  at  the  positive  end.  The  derived  solution-induced  electric  field 
strength  normal  to  the  surface  at  the  negative  end  exceeds  1 .0  V/nm,  corresponding  to  a 
total  field  strength  of  about  1 .6  V/nm  when  including  the  external  field.  This  value  is 
again  in  excellent  agreement  with  the  derived  values  for  the  concentrated  solution  (see 
value  in  Table  1,  obtained  with  the  use  of  Eq,  (6)).  The  field  is  slightly  lower  at  the 
positive  end.  The  large  surface  fields  at  the  polar  ends  of  the  drop  are  not  solely  a 
reflection  of  the  ionic  charge  density,  but  are  attributable  to  both  aligned  formamide 
dipoles  as  well  as  field-induced  polarization  and  separation  of  the  ionic  charges.  Similar 
surface  fields  are  deduced  for  the  pure  formamide  droplet  (see  Fig.  5)  where  the  absent 
contribution  to  the  field  coming  (in  the  case  of  salt  solution)  from  the  solvated  ions  is 
compensated  by  an  increased  field  due  to  a  higher  formamide  dipole-field  alignment.  As 
aforementioned,  the  above  indicates  that  in  solution  drops  the  presence  of  solvated  ions 
frustrates  the  alignment  of  the  dipolar  solvent  molecules. 

The  observed  preference  of  negative  ions  at  the  interface  under  the  influence  of  an 
external  electric  field  is  not  as  pronounced  as  the  corresponding  findings  at  aequous 
interfaces  under  field- free  conditions,  where  large  polarizable  anions  were  particularly 
preferred  over  cations  in  the  interface  layer.37, 38  In  the  present  simulations  the  anions  (F) 
are  treated  as  a  point  charge,  and  it  is  possible  that  the  inclusion  of  a  polarization  term  in 
the  interaction  potential  may  enhance  the  anionic  presence;  indeed  the  role  of  the  anionic 
polarization  term  has  been  discussed  in  the  context  of  MD  simulations  of  aequous 
solutions.  J8’  41,  42  From  the  fact  that  the  simulated  field-free  Nal/formamide  droplet 
exhibits  no  difference  between  anionic  and  cationic  densities  in  the  interface  layer  (where 
ion  depletion  occurs  as  expected  from  a  solvation  energy  consideration),  we  conclude 
that  the  minor  preference  for  the  anions  to  reside  in  the  interfacial  layer,  found  in 
simulations  under  the  effect  of  an  external  field,  is  related  to  the  interaction  of  the  ions  in 
the  cluster  with  the  applied  field.  Because  the  solvation  number  distributions  are  broader 
for  the  larger  anions,  and  because  the  ligand  binding  energies  are  weaker  (see  Table  6),  it 
may  be  expected  that  the  propensity  of  the  formamide  molecules  coordinated  with  the  1" 
to  align  with  the  field  is  higher  than  that  of  molecules  coordinated  with  the  Na+  cations, 
resulting  in  a  stronger  attraction  towards  the  high-field  region  at  the  interface.  The  lack 
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of  a  charge  preference  in  the  field-free  NaJ/formamidc  droplet  could  also  be  attributed  to 
the  nanometer-scale  droplet  confinement  effects,  surface  curvature,  and  characteristics  of 
the  liquid  solvent  (that  is,  fomamide  in  the  present  simulations  versus  water  in  previous 
interfacial  studies).  The  chemical  identity  factor  could  be  of  principal  significance, 
particularly  in  light  of  the  fact  that  the  behavior  of  aqueous  solutions  is  often  discussed 
in  terms  of  hydrogen  bonding  between  water  molecules  in  clusters  of  large  anions, 
versus  symmetric  solvation  in  the  case  of  the  cationic  systems.35  Accordingly,  it  has  been 
shown  through  MD  simulations  of  aqueous  Na+  and  C!'  clusters,  that  Na+  is  solvated 
symmetrically  in  the  center  of  the  clusters,  while  Cl'  resides  on  the  surface  of  a  water 
cluster.83  This  was  confirmed  experimentally  for  the  larger  halide  ions  in  water  dimer 
clusters  by  Johnson  and  coworkers.84  The  present  ah  initio  optimized  structures  of 
formant ide-solvated  ions,  however,  show  no  evidence  for  hydrogen  bonding  with  both 
the  form  amide- solvated  Na+  and  T  ions  located  close  to  the  centers  of  the  corresponding 
clusters. 

6.  Summary 

In  this  paper  we  presented  laboratory  experiments  and  theoretical  investigations  aiming  at 
characterizing  and  understanding  the  atomic-scale  structure,  energetics,  stability, 
dynamics,  and  fragmentation  processes  of  complex  liquids  under  the  influence  of  high 
electric  fields.  The  present  MD  simulations  of  drops  of  Nal/formamide  solutions  placed 
in  a  uniform  field  reproduce  faithfully  the  main  salient  experimental  features  observed  in 
vacuum  electrospray  ionization  experiments  of  solutions  of  similar  concentration, 
operated  in  a  single  cone-jet  mode  with  jet  dimensions  comparable  to  the  drop  size.  At 
these  conditions,  Taylor  cones  spray  charged  particles  in  a  mixed  solvated  ion/droplct 
regime.  The  MD  simulation  of  a  droplet  subjected  to  a  uniform  electric  field  of  0.625 
V/nm  shows  the  occurrence  of  droplet  fission  and  field  ion  evaporation  processes.  The 
simulation  predicts  that  ion  emission  occurs  at  the  tip  of  the  jet  that  develops  beyond  the 
neck  connecting  to  the  Taylor  cone  (see  Figs.  1,  6d  and  6e).  This  contrasts  with 
continuum  electrohydrodynamic  theory  where  emission  from  the  neck  region  is 
predicted.22,  24  The  results  from  the  simulation  are  consistent  with  the  present 
experimental  findings  that  most  of  the  field  evaporated  ions  are  formed  with  energies 
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near,  or  lower,  than  those  of  the  charged  droplets.  The  surface  normal  field  strength 
determined  from  the  simulation  near  the  point  of  charge  emission  (for  the  concentrated 
solution  drop)  is  determined  to  be  ~1.6  V/nm,  in  excellent  agreement  with  the  field 
derived  from  the  experimental  observations  using  continuum  electrostatic  theory.  The 
size  dependence  of  the  surface  electric  field  and  charge-to-mass  ratios  of  droplets  emitted 
in  the  MD  simulation  are  fully  consistent  with  the  experimentally  derived  maximum 
surface  field  strength  and  associated  curvature.  A  transition  in  the  electrostatic  properties 
of  the  emitted  clusters  as  a  function  of  size  is  observed  at  radii  of  ~1  nm,  corresponding 
to  a  transition  between  solvated  single  ions  (for  R  <  1  nm)  and  multi-ion  droplets  (for  R  > 
lnm).  This  crossover  is  attributed  to  the  stronger  intermolecuiar  forces  and  the 
consequently  higher  sustainable  charge  densities  in  the  solvated  single  ions. 

The  simulated  and  measured  solvation  number  distributions  for  Na+(HCONH2)n 
and  T(HCONH2)n  are  in  good  agreement.  While  the  Na+  cluster-ion  distributions  are 
governed  primarily  by  thermal  stability  and  are  independent  of  concentration,  the  F 
cluster-ion  distributions  broaden  with  dilution.  Similar  behaviour  is  observed 
experimentally  with  Kl/formamidc  solutions.  Larger  cluster  ions  that  are  thermally 
unstable  at  300  K  are  observed  and  attributed  to  evaporatively  cooled  species.  Complex 
ions  including  ion  pairs  are  observed  both  experimentally  and  in  the  simulations.  While 
multiply  charged  ions  are  not  observed  experimentally,  minor  quantities  of  multiply 
charged  ions  survive  the  -1 5  ps  flight  time  of  the  simulation  space. 

The  present  MD  simulations  reveal  that  for  a  Nal/formamide  droplet  at  field-free 
conditions,  ion  densities  are  depleted  in  the  interfacial  layer  and  that  there  is  no 
preference  of  anions  over  cations.  Upon  application  of  a  field,  ion  densities  in  the 
interfacial  layer  increase  at  locations  where  the  surface  field  strength  is  high.  A  slight 
preference  for  anions  is  observed  at  the  negative  droplet  end  following  field  deformation. 
This  is  attributed  to  weaker  solvent  molecule  bonding  to  F  in  comparison  to  Na\  as 
determined  by  the  present  ab  initio  calculations,  and  the  resulting  higher  degree  of 
alignment  of  the  formamide  dipole  in  the  vicinity  of  the  anion. 
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Table  1.  Electrospray  source  parameters  of  the  experimentally  investigated 
Nal/formamide  solutions,  /(tot)  and  /(Faraday  Cup)  correspond  to  the  total 
emission  and  Faraday  Cup  currents,  respectively,  m  (QCM)  (ng  s'1)  is  the  mass 
flow  as  measured  on  the  QCM.  Emax  and  Rmai  correspond  to  the  maximum 
electric  field  strength  and  the  associated  radius  of  curvature  at  the  jet  neck  given 
by  Eqs.  (6)  and  (7), 


Nal  mole  ratio 

Conductivity 

1:8.4 

1.64  S/m 

1:36.9 

1.26  S/m 

Extraction  Voltage  (V) 

2,000 

2,000 

2,400 

2,400 

Polarity 

+ 

- 

+ 

- 

I(total)  (nA) 

250 

464 

120 

200 

I(Faraday  Cup)  (nA) 

28 

58 

63 

60 

m  (QCM)  (ng  s'1) 

2.4 

3.0 

5.4 

4.0 

m/q  (Daltons) 

8,410 

4,930 

8,330 

6,360 

Emax  (V/nm) 

1.70 

1.64 

1.42 

1.49 

Rmax  (nm) 

2.30 

2.48 

3.29 

2.98 
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Table.  2.  Mass  distributions  for  Na+(HCONH2)n  and  Na+(NaI)(HCONH2)m  ions  sprayed 
from  concentrated  and  dilute  (1:8.4  and  1:36.9  Nal-formamide  mole  ratios, 
respectively)  N  a  1/form  amide  solutions.  1  is  the  recorded  normalized  mass 
spectral  intensity. 


n 

/ 

Concentrated 

Dilute 

m 

/ 

Concentrated 

i 

0.0018 

0.0007 

0 

0.0012 

2 

0.0139 

0.0111 

1 

0.0033 

3 

0.1437 

0.0967 

2 

0.0016 

4 

0.5467 

0.4674 

3 

0.0024 

5 

0.1751 

0.1505 

4 

0.0078 

6 

0.0479 

0.0699 

5 

0.0155 

7 

0.0195 

0.0418 

6 

0.0030 

8 

0.0090 

0.0335 

7 

0.0012 

9 

0.0024 

0.0240 

8 

0.0012 

10 

0.0009 

0.0169 

9 

0.0014 

11 

0.0002 

0.0148 

12 

0.0004 

0.0128 

13 

0.0148 

14 

0.0158 

15 

0.0123 

16 

0.0169 
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Table.  3.  Mass  distributions  for  I‘(HCONH2)n  and  r(NaI)(HCONH2)m  ions  sprayed  from 
concentrated  and  dilute  (1:8.4  and  1:36.9  Nal-formamide  mole  ratios,  respectively) 
Nal/fonnamide  solutions.  /  is  the  recorded  normalized  mass  spectral  intensity. 


I  / 


n 

Concentrated 

Diluted 

m 

Concentrated 

0 

0.0281 

0.0350 

0 

0.0073 

1 

0.2151 

0.0817 

1 

0.0141 

2 

0.2949 

0.0682 

2 

0.0460 

3 

0.1855 

0.1051 

3 

0.0457 

4 

0.1011 

0.1276 

4 

0.0221 

5 

0.0402 

0.1547 

6 

0.0846 

7 

0.0570 

8 

0.0692 

9 

0.0492 

10 

0.0399 

11 

0.0501 

12 

0.0292 

13 

0.0317 

14 

0.0169 
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Table.  4.  Mass  distributions  for  K+(HCONib)n  and  K*(KI)(HCONH2)m  ions  sprayed 
from  concentrated  and  dilute  (1:41  and  1:83  Kl-formamide  mole  ratios, 
respectively)  Kl/forniamide  solutions.  I  is  the  recorded  normalized  mass  spectral 
intensity.  Trace  amounts  of  unidentified  peaks  are  not  included. 


I  I 


n 

1/4 

volumetric 

diluted 

1/8  volumetric 
diluted 

m 

1/4  volumetric 
diluted 

1/8  volumetric 
diluted 

0 

0.0115 

0.0050 

0 

0.01 11 

0.0277 

1 

0.0031 

0.0009 

1 

0.0107 

0.0219 

2 

0.0085 

0.0060 

2 

0.0051 

0.0104 

3 

0.0106 

0.0421 

3 

0.0018 

0.0037 

4 

0.0861 

0.1745 

4 

0.0010 

0.0017 

5 

0.1510 

0.1758 

5 

0.0010 

0.0018 

6 

0.1368 

0.0827 

6 

0.0008 

0.0016 

7 

0.1091 

0.0576 

7 

0.0010 

8 

0.0959 

0.0476 

9 

0.0684 

0.0385 

10 

0.0497 

0.0330 

1 1 

0.0446 

0.0415 

12 

0.0462 

0.0377 

13 

0.0416 

0.0429 

14 

0.0336 

0.0374 

15 

0.0248 

0.0449 

16 

0.0334 

0.0405 
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Table.  5.  Mass  distributions  for  f'(HCONH2)n  and  I'(KI)(HCONH2)m  ions  sprayed  from 
concentrated  and  dilute  (1:41  and  1:83  Kl-formamide  mole  ratios,  respectively) 
Kl/formamide  solutions.  /  is  the  recorded  normalized  mass  spectral  intensity.  Trace 
amounts  of  unidentified  peaks  are  not  included. 


n 

1/4  diluted 

/ 

1/8  diluted 

m 

1/4  diluted 

I 

1/8  diluted 

0 

0.0076 

0,0059 

0 

0.0018 

0.0004 

1 

0.0095 

0.0142 

! 

0.0029 

0.0003 

2 

0.0444 

0.0664 

2 

0.0057 

0.0117 

3 

0.1608 

0.0772 

3 

0.0183 

0.0280 

4 

0.2465 

0.0839 

4 

0.0236 

0.0309 

5 

0.2032 

0.1305 

5 

0.0164 

0.0298 

6 

0.0453 

0.1072 

6 

0.0067 

0.0367 

7 

0.0433 

0.0569 

7 

0.0241 

8 

0.0392 

0.0616 

8 

0.0152 

9 

0.0253 

0.0416 

9 

0.0176 

10 

0.0206 

0.0303 

10 

0.0157 

11 

0.0294 

0.0383 

12 

0.0263 

0.0198 

13 

0.0228 

14 

0.0185 

44 


Table.  6.  Thermochemistry  of  Na+(NH2CHO)n,  n  =  1-6,  K+(NH2CHO)„,  n  =  l-6,and 
I'(NH2CHO)„,  n  =  1-5  solvated  ions  determined  with  the  indicated  ab  initio  methods. 
The  internal  energy,  Ein,t  signifies  the  298.15  K  thermal  energy  above  the  zero-point 
energy,  E:pe.  The  binding  energy,  D(MFn.i-F),  refers  to  the  loss  of  a  single  solvent 
molecule,  F=HCONH2.  Also  shown  are  thermochemical  data  at  310  K  detennined  for 
the  “classical”  clusters  bound  by  theMD  simulation  inter-atomic  potentials.  Ail  values 
are  given  in  kcal/mol. 


Quantum  Chemistry 

E:pe  Ein,  D(MFn.,-F) 

MI) 

E/nt  D(MF„.  j— 

F) 

Na+(NH2CHO)„ 

B3LYP/6-31G(d) 

1 

30.1 

3.33 

39.7 

1.84 

21.6 

2 

59.9 

6.7! 

33.3 

3.69 

19.8 

3 

89.4 

10.3 

24.1 

5.53 

16.8 

4 

118.4 

14.0 

17.2 

7.40 

14.3 

5 

149.5 

16.3 

15.7 

9.24 

11.3 

6 

178.2 

20.2 

7.2 

11.1 

9.57 

K+(NH2CHO)„ 

B3LYP/6-31G(d) 

1 

29.8 

3.46 

28.7 

2 

59.3 

7.02 

24.0 

3 

88.7 

10.6 

19.5 

4 

117.7 

14.4 

15.6 

5 

148.8 

16.8 

34.4 

6 

178.3 

20.3 

10.8 

I(NH2CHO)n 

MP2/6-3 1 1  G+(d,p)/6-3 1 1+G(df> 

1 

27.8 

3.6 

17.4 

1.85 

13.1 

2 

56.0 

7.0 

14.9 

3.70 

12.9 

3 

83.7 

10.8 

12.1 

5.54 

11.3 

4 

111. 8 

14.3 

10,9 

7.39 

10.2 

5 

140.0 

17.1 

8.4 

9.24 

8.6 

45 


Table.  7.  Survival  fractions  for  selected  Na+(NH2CHO)„  ions  evaporated  with 
characteristic  temperatures  of  300,  350  and  400  K.  The  survival  fractions  are  listed  for 
infinite  observation  times,  t,  flight  times  through  mass  spectrometric  experiment  (MS), 
and  the  MD  simulation  period  (MD). 


/ 

n 

4 

5 

6 

T=  300  K 

OO 

0.927 

0.503 

8.0  x  10‘s 

MS  (26, 28,  31  ps) 

0.994 

0.662 

1.0  x  104 

MD  ( 1 5  ps) 

T  =  350  K 

1.000 

0.903 

7.2  x  10'3 

co 

0.702 

0.148 

0,000 

MS  (26,  28,  31  ps) 

0.920 

0.279 

0.000 

MD  ( 1 5  ps) 

T  =  400  K 

1. 000 

0.601 

8.0  x  10’4 

oo 

0.332 

0.035 

0.000 

MS  (26,  28.  31  ps) 

0.718 

0.070 

0.000 

MD(15  ps) 

1.0 

0.273 

0.000 
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Table.  8.  Survival  fractions  for  selected  I'(HCONH2)„  ions  evaporated  with  characteristic 


temperatures  of  300,  350  and  400  K.  The  survival  fractions  are  listed  for  infinite 

observation  times,  f,  flight  times  through  mass  spectrometric  experiment  (MS),  and  the 

MD  simulation  period  (MD), 

/ 

n 

3 

4 

5 

00 

T=  300  K 

0.859 

0.295 

0.036 

MS  (29,  32,  34 
us) 

0.957 

0.548 

0.091 

MD  (15  ps) 

1,000 

0.974 

0.407 

00 

T  =  350  K 

0.620 

0.088 

0.006 

MS  (29,  32,  34 

MS) 

MD(15  ps) 

0.811 

0.224 

0.016 

1.000 

0.837 

0.121 

00 

T =400  K 

0.341 

0.023 

0.000 

MS  (29,  32,  34 
ps) 

MD  (15  ps) 

0.570 

0.067 

0,001 

1.0 

0.560 

0.024 
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Figure  Caption: 

Fig.  1 .  Schematic  representation  of  an  electrospray  cone-jet. 

Fig.  2.  Schematic  diagram  of  the  mass  spectrometric  experiment. 

Fig.  3.  Sample  positive  and  negative  ion  mass  spectra  of  the  concentrated  (28.4%  wt.) 
Nal-formamide  solution,  Solvated  ion  solvation  number  assignments  are  shown. 

Fig.  4.  Retarding  potential  measurements  of  solvated  ion  and  droplet  (corresponding 
to  high  pass,  mtq  >  400)  energy  distributions,  The  capillary  potential  was  biased  at 
550  V,  and  the  analyzer  was  at  ground. 

Fig.  5.  MD  simulation  of  a  1 0  nm  fonnamide  droplet  prior  to  and  after  the  application 
of  an  electric  field.  When  a  critical  field  is  exceeded,  the  droplet  lengthens  to  another 
stable  configuration.  The  calculated  critical  field  for  a  10  nm  droplet  assumes  a  surface 
tension  of  0.046  N/m.  The  aspect  ratios,  r|,  for  the  various  configurations  are 
specified  in  the  figure. 

Fig.  6.  MD  simulation  snapshots  of  a  10  nm  diameter  droplet  of  a  1:16  ratio 
Nal/Formamide  solution  during  different  stages  of  elongation  due  to  the  application  of 
a  uniform  external  electric  field  of  0.625  V/nm.  Na4  and  I"  are  shown  as  red  and 
yellow  spheres  respectively.  Time  with  respect  to  field  application;  a)  0  ps;  b)  400  ps; 
c)  700  ps;  d)  1 ,1 00  ps;  e)  !  ,1 25  ps.  The  aspect  ratios,  rp  for  the  various  configurations 
are  specified  in  the  figure. 


Fig.  7.  Properties  of  the  charged  clusters  emitted  from  the  droplet  described  in  Fig.  5 
as  a  function  of  cluster  radius:  (a)  the  estimated  electric  field  at  the  droplet  surface,  (b) 
the  charge  to  mass  ratios  and  (c)  the  observed  (MD)  charge  divided  by  the  predicted 
charge  for  Rayleigh  instability  and  fission.  The  added  dashed  and  thin  solid  lines 
compare  MD  results  to  estimates  based  on  the  experimental  electrospray  results  (from 
Table  I ,  see  text  for  details). 
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Fig.  8.  Number  densities  and  induced  electric  field  of  the  1:16  ratio  Nal/Formamide 
droplet  along  the  field  direction  for  (a)  no  external  field,  (b)  elongating  in  constant 
external  0.625  V/nm  field  before  emitting  ions  and  (c)  during  the  time  of  initial 
cluster-ion  emissions.  The  profiles  correspond  to  the  respective  snapshots  in  Fig.  5a-c. 
The  droplet  components  are  designated  as  in  (a)  while  the  relative  densities  of  Na  to  I 
and  the  electric  field  are  as  given  in  (b)  and  (c).  The  internal  field  of  the  droplet 
(generated  by  induced  surface  charges)  is  shown  as  lighter  lines  and  fluctuates  around 
-  -0.6  V/nm.  The  aspect  ratios,  t|,  for  the  various  configurations  are  specified  in  the 
figure. 

Fig.  9  Charge  abundances  of  ion-containing  clusters  emitted  from  a  1:8  ratio 
Nal/Formamide  droplet  (29%  by  weight  of  dissolved  Nal)  placed  in  a  uniform  external 
field.  Results  are  shown  for  two  simulations  having  external  field  strengths  of  1 .0  and 
1 .5  V/nm.  The  observation  occurred  at  ±200  nm  from  the  droplet. 

Fig.  10,  Ion  composition  of  singly  charged  formamide/Nal  clusters  emitted  from  a  1 :8 
ratio  Nal/Formamide  droplet  (29%  by  weight  of  dissolved  Nal)  placed  in  a  uniform 
external  field  strengths  of  1.5  V/nm  (top)  and  1.0  V/nm  (bottom).  The  observation 
occurred  at  ±200  nm  from  the  droplet. 

Fig.  1 1 .  Solvation  number  distribution  of  single-ion  clusters  emitted  from  a  1 :8  ratio 
Nal/Formamide  droplet  (29%  by  weight  of  dissolved  Nal  )  placed  in  a  uniform  external 
field  of  1 .5  V/nm.  The  observation  occurred  at  ±200  nm  from  the  droplet.  MD  results 
are  compared  with  experiment.  MD  results  corrected  for  30  ps  flight  time  (adjusted 
MD)  are  also  shown  (see  text  for  details). 
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relative  abundance 
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Fig  4. 
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field  induced  elongation  offormamide  droplet 
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Fig  5. 
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Fig  7. 
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